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SUMMARY 


A  topic  of  ongoing  Interest  to  the  Air  Force  Is  the  response  of 
shell  structures  to  thermal  loads.  These  thermal  loads  may  be  caused  by 
aerodynamic  heating  on  aircraft  and  missiles,  by  laser  weapons,  or  by 
radiation  effects  on  space  structures.  In  all  cases.  It  Is  desirable  to 
determine  the  thermomechanical  response  of  shell-llke  bodies  that  model 
various  components  of  aircraft  structures. 

To  determine  the  complete  thermomechanical  response  of  a  shell-llke 
body  to  thermal  loads.  It  Is  necessary  to  calculate  both  the  thermal  and 
mechanical  fields.  The  problems  with  analytically  solving  the  three- 
dimensional  thermal  and  mechanical  equations  for  general  structural 
geometry  and  general  loading  are  currently  Insurmountable.  Therefore, 
alternative  methods  of  solution  are  required.  One  alternative  Is  to 
obtain  a  numerical  solution  using  large  finite  element  codes.  Although 
obtaining  a  numerical  solution  of  this  kind  is  within  current  technol¬ 
ogical  capability,  there  are  two  Important  disadvantages  with  this 
method.  First,  It  may  be  necessary  to  use  many  elements  and  many  time 
steps  to  calculate  the  solution  of  a  dynamic  problem.  This  means  that 
use  of  the  method  as  a  design  tool  when  a  number  of  parameters  must  be 
varied  would  be  prohibitively  expensive.  Second,  the  method  necessarily 
calculates  many  details  of  the  thermal  and  mechanical  fields  that  are 
not  of  particular  Interest.  For  example,  for  shells  it  Is  sufficient  to 
obtain  limited  Information  such  as  the  resultant  force  and  moment  applied 
to  the  shell.  Consequently,  the  details  of  the  stress  distribution 
through  the  thickness  of  the  shell  are  not  needed. 

In  this  research,  we  used  an  alternative  method  of  solution  that 
judiciously  models  the  structure  by  an  appropriate  shell  theory  having 
equations  that  are  considerably  simpler  than  those  of  the  three- 
dimensional  theory.  Such  a  theory  focuses  attention  solely  on 
quantities  of  primary  Interest.  Consequently,  even  If  It  Is  necessary 


to  develop  a  numerical  solution  of  the  shell  equations,  the  computer 
time  Is  efficiently  used,  the  results  are  easily  Interpreted,  and  the 
method  can  be  used  for  design  purposes. 

The  objective  of  this  research  was  to  analyze  the  thermomechanical 
response  of  a  conical  shell  (Figure  1)  that  Is  a  model  for  aircraft  and 
missile  nose  cones.  Because  the  conical  shell  has  a  converging 
geometry,  the  shell  near  Its  tip  Is  necessarily  "thick'*  even  though  the 
shell  near  Its  base  may  be  "thin.”  Therefore,  It  Is  Important  to  model 
the  conical  shell  with  a  theory  that  accurately  Incorporates  the 
geometrical  details  of  this  crucial  tip  region.  We  used  recent  devel¬ 
opments  In  the  Cosserat  theory  for  the  thermomechanical  response  of 
shells.  We  based  our  development  on  the  Cosserat  theory  mainly  because 
this  theory  accurately  models  the  geometry  of  the  shell  and  Is 
sufficiently  general  enough  to  Include  the  Important  effects  of  the 
steep  temperature  profile  through  the  thickness  of  the  shell  (e.g., 
tangential  shear  deformation,*  higher  order  displacement  and  temperature 
profiles  through  the  thickness  of  the  shell)  without  the  complexity  of 
the  complete  three-dimensional  theory. 

This  research  t«s  divided  Into  four  parts,  each  of  which  has  been 
written  as  a  technical  paper  that  has  been  submitted  for  publication.  A 
copy  of  each  paper  Is  Included  as  an  appendix: 

Appendix  A:  A  Uniqueness  Theorem  for  Thermoelastic  Shells 
with  Generalized  Boundary  Conditions 

Appendix  B:  A  Nonlinear  Constrained  Theory  of  Shells  that 
Includes  Tangential  Shear  Deformation 

Appendix  C:  On  the  Determination  of  Certain  Constitutive 
Coefficients  for  Theinnoelastlc  Shells 

Appendix  D:  Heat  Conduction  In  Flates  and  Shells  with 
Emphasis  on  a  Conical  Shell. 

*We  use  the  term  tangential  shear  deformation  Instead  of  the  usual  term 
"transverse'*  shear  deformation  because  It  Is  more  descriptive  when 
considering  nonlinear  deformation  of  a  shell  as  opposed  to  linear 
deformation  of  a  plate. 
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FIGURE  1 


CONICAL  SHELL  OF  CONSTANT  THICKNESS  h, 
TIP  RADIUS  R, .  AND  BASE  RADIUS  Rj 
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Basically,  Appendices  A-C  deal  with  certain  theoretical  Issues  that 
were  observed  while  developing  the  general  equations  to  analyze  the 
therraomechanlcal  response  of  shells,  whereas  Appendix  D  deals  with  the 
main  Issue  of  the  conical  shell.  Appendix  A  entails  a  generalization  of 
the  uniqueness  theorem  for  linear  thermoelastic  shells  to  Include 
mechanical  contact  with  a  linear  elastic  medium  and  thermal  convection 
on  both  the  major  surfaces  of  the  shell  and  Its  boundary  curve. 

Appendix  B  Is  concerned  with  developing  a  nonlinear  constrained  theory 
of  shells  that  Includes  tangential  shear  deformation.  This  topic  was 
considered  because  tangential  shear  deformation  Is  expected  to  be 
Important  when  a  shell  Is  loaded  with  a  steep  through- the- thickness 
temperature  gradient.  Appendix  C  Is  concerned  with  determining  certain 
thermal  and  mechanical  constitutive  coefficients  for  thermoelastic 
shells  by  direct  comparison  with  exact  solutions.  It  was  observed  that 
better  comparison  with  simple  exact  solutions  could  be  obtained  by 
specifying  values  for  certain  coefficients  that  are  different  from  those 
previously  proposed. 

In  addition  to  considering  general  aspects  of  thermomechanical 
response  of  shells.  In  Appendix  D,  we  considered  heat  conduction  In 
rigid  shells  with  particular  emphasis  on  a  conical  shell.  Specifically, 
we  used  shell  equations  based  on  the  theory  of  a  Cosserat  surface  to 
determine  the  average  (through  the  thickness)  temperature  and  tempera¬ 
ture  gradient  In  rigid  shells.  Attention  was  focused  on  rigid  shells 
because  when  the  strain  rates  In  a  deformable  shell  are  small  enough, 
the  thermal  and  mechanical  problems  uncouple  In  the  sense  that  the 
temperature  fields  can  be  determined  from  equations  for  heat  conduction 
In  rigid  shells.  Once  these  temperature  fields  are  known,  they  can  be 
used  together  with  constitutive  equations  to  calculate  thermal  loads 
which  cause  mechanical  deformation. 

To  develop  confidence  In  the  Cosserat  theory  for  both  the  thin- 
shell  limit  (which  models  the  base  of  the  conical  shell)  and  the  thick- 
shell  limit  (which  models  the  tip  of  the  conical  shell),  we  considered  a 
number  of  problems  for  plates,  circular  cylindrical  shells,  spherical 


shells,  and  finally  a  conical  shell  (Figure  1).  It  was  shown  that  by 
appropriately  modifying  the  constitutive  equations.  It  Is  possible  to 
Include  enough  geometrical  features  of  the  shell  to  predict  relatively 
accurate  results  even  In  the  thick-shell  limit. 

Finally,  a  specific  problem  of  a  conical  shell  of  uniform  thickness 
was  considered.  The  heat  flux  on  the  outer  surface  of  the  shell  was 
taken  to  be  constant,  all  other  surfaces  of  the  shell  were  Insulated, 
and  the  shell  was  Initially  at  uniform  temperature.  Using  the  Cosserat 
equations,  we  were  able  to  obtain  an  analytical  solution  for  the  time- 
dependent  average  temperature  and  temperature  gradient.  This  solution 
predicted  that  the  average  temperature  gradient  was  quite  severe  near 
the  tip  of  the  shell  and  developed  its  maximum  value  within  a  short 
time.  Since  the  average  temperature  gradient  produces  a  thermal  moment 
on  the  shell,  we  expect  the  bending  to  be  most  severe  near  the  shell's 
tip. 

Future  research  should  concentrate  on  determining  the  deformation 
In  the  conical  shell  produced  by  these  thermal  fields. 
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A  UNIQUENESS  THEOREM  FOR  THERMOELASTIC  SHELLS 
WITH  GENERALIZED  BOUNDARY  CONDITIONS 


(Submitted  for  publication  to  Quarterly  of 
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ABSTRACT 


The  uniqueness  of  the  solution  of  Initial,  mixed  boundary  value 
problems  for  linear  thermoelastic  sheila  Is  reconsidered  within  the 
context  of  recent  developments  In  the  thermomechanical  theory  of  a 
Cosserat  surface  [4].  Fairly  general  boundary  conditions  are  considered 
that  allow  mechanical  contact  with  linear  elastic  media  and  thermal 
radiation  on  the  boundary  curve  of  the  Cosserat  surface  and  on  the  major 
surfaces  of  the  shell. 
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Introduction 


Motivated  by  the  application  of  shell  theory  to  contact  problems, 
we  generalize  the  usual  boundary  conditions  to  Include  mechanical 
contact  with  linear  elastic  media  and  thermal  radiation.  These 
generalized  mechanical  and  thermal  boundary  conditions  are  applied  on 
the  boundary  curve  of  the  shell  as  well  as  on  Its  major  surfaces.*  In 
this  paper,  the  uniqueness  proof  Is  reconsidered  within  the  context  of 
these  general  boundary  conditions. 

We  recall  that  the  theory  of  a  Cosserat  surface  [1]  has  been  well 
established  as  a  particularly  useful  model  of  a  shell-llke  body  that, 
broadly  speaking.  Is  a  three-dimensional  body  that  Is  "thin”  In  one  of 
Its  dimensions.  Recent  developments  In  the  theory  of  continuum  thermo¬ 
dynamics  [2,3]  have  provided  the  theoretical  framework  for  developing  a 
general  thermomechanical  theory  of  a  Cosserat  surface  [4].  Such  a 
theory  admits  finite  numbers  of  directors  d  (N  -  1,2...)  and  tempera- 
ture  field  6^  to  provide  limited  Information  about  the  variation  through 
the  thickness  of  the  shell  of  the  deformation  and  temperature  field, 
respectively. 

For  less  general  boundary  conditions  than  those  considered  here,  a 
uniqueness  theorem  has  been  proved  for  the  linear  Isothermal  theory  of 
shells  [1,  Sec.  26]  and  for  small  motions  superposed  on  a  large  deforma¬ 
tion  within  the  context  of  a  thermoelastic  theory  of  shells  that  admits 
a  single  temperature  field  [5] .  A  uniqueness  theorem  for  a  linear 
thermoelastic  theory  of  shells  that  admits  two  temperature  fields  and 
two  energy  equations  has  also  been  proved  [6],  Including  radiation  on 
the  major  surfaces,  but  not  on  the  boundary  curve.  None  of  these 
uniqueness  proofs  considers  mechanical  contact  with  elastic  media. 


Recall  [1]  that  the  three-dimensional  boundary  conditions  applied  to 
the  major  surfaces  of  the  shell  are  Incorporated  Into  the  field 
equations  and  therefore  are  not  considered  to  be  boundary  conditions 
in  shell  theory. 


A-3 


Although  the  structure  of  this  latter  theory  Is  different  from  that 
most  recently  developed  for  Cosserat  surfaces  [A],  the  linearized 
equations  can  be  placed  in  a  one-to-one  correspondence  with  those  of  a 
theory  that  admits  a  single  director  d  and  two  temperature  fields  9,  4>* 
It  therefore  follows  that  the  previous  uniqueness  theorem  [6]  applies  to 
solutions  within  the  context  of  the  new  linearized  theory  [A]. 

In  the  following  sections,  we  state  sufficient  conditions  to  prove 
uniqueness  of  the  solution  of  the  equations  of  the  coupled  thermoelastic 
theory*  that  admits  the  generalized  boundary  conditions  and  that 
considers  Inhomogeneous,  anisotropic  elastic  shells.  Specifically, 
Section  2  records  the  basic  equations  of  the  linear  theory  of  a  Cosserat 
surface,  and  Section  3  discusses  the  generalized  boundary  conditions. 
Finally,  in  Section  A  we  state  and  prove  the  uniqueness  theorem. 

2.  Basic  Equations 

Let  the  material  points  of  the  Cosserat  surface  C  be  identified  by 
means  of  a  system  of  convected  coordinates  6°^  (a  ~  1,2)  and  let  the  two- 
dimensional  region  occupied  by  the  material  surface  in  the  present 
configuration  at  time  t  be  denoted  by  c.  Further,  let  the  vector  valued 
function  r  define  the  position  of  a  material  point  of  the  surface  C  and 
at  each  such  point  define  the  vector  valued  function  d,  called  the 
director,  and  the  two  temperature  fields  9  and  ^  each  referred  to  the 
present  configuration.  Then,  a  thermomechanical  process  of  the  Cosserat 
surface  is  defined  by 

r  -  rCO*.!;)  ,  d  -  d(9‘.t)  .  [a,  a,  d]  >  0  (2.1a, b.c) 

9  -  9(9®, t)  ,  (9  >  0)  ,  ♦  -  ♦(9®,t)  ,  (2.1d,e,f) 


The  uniqueness  theorem,  with  its  modification  for  the  generalized 
boundary  conditions,  is  clearly  applicable  to  the  purely  mechanical 
theory  as  well  as  the  purely  thermal  theory. 


A-A 


where  the  tangent  vectors  a  and  unit  normal  vector  a.  are  defined  by 

~o  ~3 


a  -  -  ,  a  •  a.  ■  0  ,  a,  •  a,  -  1  ,  a 

~a  *  ~o  ~3  *  ~3  ~3 


fSl  £2  ^3^  >  ° 
(2.2a,b,c,d) 


and  the  condition  (2.1c)  ensures  that  the  director  Is  nowhere  tangent  to 
c.  The  condition  (2. 2d)  ensures  that  a  (1  ■  1,2,3)  are  linearly  inde- 
pendent  vectors  forming  a  right-handed  coordinate  system.  Thus,  we  may 
Introduce  a  set  of  reciprocal  vectors  a^  such  that 


j  .j 
£1  *  £  " 


(2.3) 


where  6;,  is  the  Kronecker  symbol.  The  velocity  v,  director  velocity  w, 

1  ^  -s. 

and  temperature  gradients  g  and  g,  may  now  be  defined  by 


V  -  r  , 


£  -  £  •  Jl  -  £  »  Jii  “  ♦,£ 


(2.4a,b,c,d) 


where  a  superposed  dot  denotes  time  differentiation  holding  0*^  fixed  and 
a  comma  denotes  partial  differentiation  with  respect  to  the  coordinates 
9°^.  In  the  reference  configuration,  we  assume  that  the  shell  Is  at 
uniform  temperature  ®.  Then,  the  reference  values  of  the  various 
kinematic  quantities  may  be  denoted  by 


r-R  ,  d-D  ,  a,  -  A, 


1/2  1/2 

a  “A 


(2.5a,b,c,d) 


,  0  -  0  , 


(2.5e,f) 


1/2  a 

where  R,  D,  ^  and  A  depend  on  the  coordinates  0  only. 


Throughout  this  text,  we  use  the  usual  summation  convention  over 
repeated  Indices.  Greek  Indices  have  the  range  (1,2)  and  Latin  Indices 
have  the  range  (1,2,3). 


.> 


Let  P,  bounded  by  the  closed  curve  dP,  denote  the  region  occupied 
by  an  arbitrary  material  portion  of  the  surface  c  In  the  present 
configuration  and  let  v  be  the  unit  outward  normal  to  dP.  Using  the 
notation  of  [7]  and  referring  to  the  present  configuration,  we  define 
the  following  quantities:  the  positive  mass  density  (mass  per  unit  area 
of  P)  p  -  p(6*,t);  the  contact  force  n  ■  n(0®,t;  v)  and  the  contact 
director  couple  m  ■  m(0*,t:  v),  each  per  unit  length  of  the  curve  SP; 
the  specific  (per  unit  mass  of  P)  assigned  force  f  >  f(0°‘,t)  and 
specific  assigned  director  couple  jt  >  1(0^, t);  the  Intrinsic  director 
couple  k  ■  k(0*,t)  per  unit  area  of  P;  the  Inertia  coefficients  y^  ■ 
y^(0“)  and  y^  -  y2(0«)  which  are  Independent  of  time;  the  specific 
entropies  q  ■  Ti(0“,t)  and  -  Tii(0*»t);  the  specific  Internal  rates  of 

production  of  entropy  5  -  5(9®»t),  -  5i(0®,t)  and  -  5^(0*. t);  the 

entropy  fluxes  k  ■  k(0®,t;  v)  at^  k,  ■  ki(0®,t;  v)  each  per  unit  length 
acting  across  the  curve  dP;  the  specific  external  rates  of  supply  of 
entropy  s  ■  a(0®,t)  and  s^  “  s^(0®,t);  the  specific  Internal  energy 
e  ■  e(0®,t);  and  the  specific  Helmholtz  free  energy  (|>  »  (|>(0“,t)  ■  e  -  0ti 
- 

For  the  linear  theory,  It  Is  convenient  to  Introduce  the  displace¬ 
ment  vector  u  and  director  displacement  6  by  the  equations 

r*R+u  ,  d»D  +  6  •  (2.6a,b) 

^  ^  ^  M  ^ 

If,  In  the  reference  configuration,  the  shell  Is  free  of  assigned  fields 
and  contact  forces  and  director  couples,  then  for  the  linear  theory  we 
assume  that  In  the  present  configuration,  the  displacements  u,  5  and  the 
temperatures  (0  -  ®)  and  ^  are  of  order^  £(0  <  e  «  1)  and  that 
quadratic  terms  In  these  quantities  may  be  neglected  relative  to  linear 

*Thls  notation  differs  from  that  used  in  [1,4]. 

^The  temporary  use  of  this  symbol  for  the  small  parameter  should  not  be 
confused  with  the  use  of  the  same  symbol  elsewhere  for  the  Internal 
energy . 
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wl' 


terms.  With  this  background,  we  record  [1]  expressions  for  the 


kinematic  quantities  A  B  A,  ,  and  for  the  strains  e  In 

ap  ap  la  ap  i  la 

the  forms 


A  _  “  A  •  A_  ,  B  ■  A  •  A-  ,  A,  “  A  •  D 
op  ~a  ap  -^a.P  ~3  la  ~1  ~,a 


(2.7a,b,c) 


‘Ua  +  A-'U  ) 

ap  2  ~a .  ~,p  ~»a 


(2.7d) 


Y  *A  ‘S+D^u  ,  Y,  *A_»5-(D«  A*)(A,  •  u  )  ,  (2.7e,f) 

'a  ~a  ~  ~  ~,a  *  '3  ~3  ~  ~  -''•/s.3  ~,a 


ic„“A  •6o+Da*u  >  <■»„  ■  A-  •  6  -  (A  •  D  ^)(A,  •  u  ) 

ap  ~a  ~,p  ~,p  ~,a  3a  ~3  ~,a  ~  ~,a  ~3  ~,a 


(2.7g,h) 


where  for  the  linear  theory  all  tensor  quantities  are  referred  to  the 


base  vectors  A. . 

~1 


With  suitable  continuity  assumptions,  It  can  be  shown  that  [1,4] 


a  a 

n  -  N  V  ,  m  -  M  V 

(X  ^  ^  (X 


(2.8a,b) 


k  .  £  .  V  -  p“  •  V  -  p“ 


(2.8c,d) 


where  v  ■  v  «  A  are  the  components  of  the  normal  vector  v  and  where 
o  ~  ~a  ~ 

N  ,  M*,  £,  and  £^  are  Independent  of  v.  Further,  with  reference  to  the 


energy  equation,  the  specific  external  rates  of  heat  supply  r  and  r^; 


and  the  heat  flux  vectors  q  and  q,  are  defined  by 


r  -  es  ,  -  ((»s^  ,  q  *  9p  . 


(2.qa,b,c,d) 


Now  the  local  forms  of  the  basic  equations  may  be  recorded  as 


^  1/2  1/2  .  „ 
\-pa  ”P_A  >0  , 


(2.10a) 


A(v  +  w)  -  +  (A^^^  n“)^^  , 


(2.10b) 


Vi"! 


‘  ^  Cvlv  VAs‘‘o‘n 


a 


(2.10c) 


1  •  -u  2  .1/2  .  ^  ,.1/2 

Hy  V  +  y  w)  -  XI  -  A  k  +  (A 

X^  -  X(8  +  O  -  (A^^^  p“)  „  , 

»o 

Xn^  -  X(s^  +  ^l>  “  Pl\a  ' 


m“) 


(2.10d) 

(2.10e) 


where  Is  the  reference  value  of  p.  Equation  (2.10a)  represents 

conservation  of  mass,  (2.10b)  represents  the  balance  of  linear  momentum, 

(2.10c)  represents  the  balance  of  director  momentum,  and  (2.10d,e) 

d  d 

represent  balances  of  entropy*  Referring  the  quantities  N  >  M  to 
the  base  vectors  A. ,  we  may  write* 

N®  -  N^®  A,  ,  k  -  A,  .  M®  -  M^®  A,  (2.11a,b,c) 

Using  these  definitions,  the  results  of  the  balance  of  angular  momentum 
become  [1] 


N^®  -  N®^  -  N^®  -  k®  -  M®®  A^O  , 


„3«  .  o  -flt  ,  3  .  .3  no  .o  3a 

N  ■Dk-Dk+  A*a  M  -  A»a  M  , 


(2.12a) 

(2.12b) 


where  tensor  quantities  with  superscripts  are  contravarlant  or  mixed 
tensor  quantities  referred  to  A^  and  A  . 

Using  appropriate  constitutive  equations  for  an  elastic  material 
and  demanding  that  the  balance  of  energy  be  Identically  satisfied  for 
all  processes.  It  can  be  shown  [4]  that 


N®P. 


S<|> 

Po^ 


,  k 


aP 


0  M^®  -  D 

°  ’  Po  * 


^  -  ae  •  ^1  “  ^  * 


(2.13a,b,c) 

(2.13d,e) 


This  notation  Is  consistent  with  [7]  but  differs  from  that  used  In 

[l.AJ. 
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and  that 


i 


+  P  •  &  +  £i  •  &i  -  0 


(2.14) 


where  for  the  linear  theory,  we  note  that  p  can  be  replaced  by  in 
expressions  of  the  type  (2.13)  and  (2.14). 

Specifically,  for  an  Inhomogeneous,  anisotropic  shell  we  assume 


2  p^  <1,  -  2  p^  4,^  -  2  Y3(e  -  )  -  2  e^p(9  -  ) 

-  2  C®^  t  -  93(9^  -  2  9  )  -  4)^  -  2  P5  9 

2  "o  *1  -  '=f 'as  'tS  'aS  't6 

*  ''a  ''S  '=f  '3.  'sp  ■"  2  'f  'as  ''3  • 

•  ■  'f  ®,S  ’  "l  •  -  *=7^  ♦.S  • 


(2.15a) 


(2.15b) 


(2.15c,d) 


■>0  ® ®.p  *  ♦,»  ♦,?  ♦  ‘2 


"o  ^1  -  p„  -  -  >>2  ♦  . 


where  the  quantities 


cf’’*CH.1.2)  .  cj;% 


(2.15f) 


1,2,. ..,7)  , 


®4  »  Pq  t  P3  »  P^  »  P3  >  b^  ,  b3  , 


are  functions  of  the  variables 


V  -  {9  ,  A^p  ,  B^p  ,  ,  }  , 


(2.16a,b) 


(2.16c) 


(2.17) 
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and  the  tensors  (2.16a,b)  satisfy  the  symmetry  restrictions 


^a3Y6  _  j,Poy6  _  ^o36y  _  ^Y6»^ 


(2.18a) 


oPy6  YfiaP 
^2  ^2 


-  c5“(N  -  1-4,  6-7) 


(2.18b,c) 


The  forms  of  the  constitutive  equations  (2.15)  are  chosen  Co  be  similar 
to  chose  used  in  Che  linearized  theory  [4],  but  with  appropriate 
generalizations  Co  allow  for  anisotropic,  inhomogeneous  response  and  to 
satisfy  the  reduced  energy  equation  (2.14)  without  approximation. 
Specifically,  equations  (2.15c-f)  satisfy  the  nonlinear  form  of  Che 
reduced  energy  equation  (2.14),  with  replaced  by  p  in  each  of  the 
expressions.  In  particular,  we  note  from  (2.15e)  that  is  of  order  e 
and  therefore  may  be  neglected  in  equation  (2.10d).  Furthermore,  we 
note  Chat  the  constitutive  equations  (2.15)  automatically  satisfy 
restrictions  for  modeling  a  shell  chat  has  S3niimecric  response  about  its 
reference  surface  [1,  Sec.  13  y1»  Now,  with  the  help  of  the  restric¬ 
tions  (2.13),  (2.18)  and  the  constitutive  equations  (2.15),  we  conclude 
Chat 

'.9  .  ^  ^.8  .  j.8(,  .  3,,  _  (2.1,., 

k®  -  c,®  Yp  ,  k^  -  a.  Y3  +  e.p  “  8.(6  “  ®)  »  (2.19b,c) 

M®»  -  Cfl'®  k.^5  -  of  ♦  .  -  of  kjp  ,  (2.l6d.,) 

Po^  ’  Po  ^3  ^3^®  -  ®)  +  P5  ,  -  cf  x^p  +  4.  , 

(2.19f,g) 

2  P^€  -  2  p^4,^  +  2  Po  Y3  +  2  cf  e^p  +  P3  9^  +  P^  (2.19h) 


Finally,  we  recall  [4]  that  for  such  an  elastic  shell  the  only 
nontrivial  statements  of  Che  second  law  of  thermodynamics  take  the  forms 


(2.20a) 


pQ  *  ^1  +  £  •  S  +  £l  *  Si  ^  °  • 

0(t)  -  9^  >  0  whenever  e(t)  -  >  0  .  (2.20b) 

The  statement  (2.20a)  corresponds  to  the  classical  heat  conduction 
Inequality  and  Is  assumed  valid  for  all  equilibrium  displacement  and 
temperature  fields.  Further,  the  statement  (2.20b)  Is  assumed  valid 
when  the  Cosserat  surface  Is  at  rest  and  the  three-dimensional  tempera¬ 
ture  field  Is  spatially  uniform  so  that  0  -  9(t),  ^  0.  In  (2.20b), 

and  0^  correspond  to  the  Internal  energy  and  uniform  temperature  of  the 
shell  during  some  period  of  time  up  to  t^  when  the  shell  has  been  at 
rest  and  In  thermal  equilibrium. 

3 .  Boundary  and  Initial  Conditions 

This  theory,  which  Is  developed  by  direct  approach,  may  be  brought 
into  a  one-to-one  correspondence  with  the  three-dimensional  theory  by 
assuming  that  the  position  vector  £  of  a  point  In  the  shell  and  the 
temperature  field  0  admit  the  representations. 

£*  •  £*(9®,9^,t)  -  r(0“,t)  +  0^  d(0®,t)  ,  (3.1a) 

9*  -  9*(9®,9^,t)  -  9(9®, t)  +  9^  (K0®,t)  ,  (3.1b) 

where  0  Is  a  coordinate  through  the  thickness  of  the  shell.  Without 
loss  In  generality,  we  may  define  the  top  surface  dP'*'  of  the  shell  by 
0  ■  h/2,  where  h  Is  a  constant  having  the  dimensions  of  length;  then  we 
may  write  the  displacement  u^  and  temperature  difference  (0'*’  -  ®)  on  5P^ 
by 


u^  -  u  +  y  6  ,  (9^  -  ®)  ■  (9  -  0)  +  Y  4)  .  (3.2a,b) 

Similarly,  we  may  define  the  bottom  surface  5P  of  the  shell  by  0  ■ 

-h/2  and  write  the  displacement  u  and  temperature  difference  (9”  -  ®) 
on  ap“  by 
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(3.3a,b) 


u  -11-^6  ,  (0  -®)-(0-@)-Y^ 


In  this  regard,  we  note  Chat  Che  specification  that  the  major  surfaces 
of  the  shell  are  defined  by  0  *  i  h/2,  where  h  is  constant,  does  not 

restrict  the  theory  to  that  associated  with  a  shell  of  uniform  thickness; 
this  is  because  Che  director  D  In  Che  reference  configuration  may  be 
specified  as  a  function  of  0^« 

With  the  help  of  (2.12)  and  the  constitutive  equations  (2.15c,d,f), 
(2.19a-g)  and  the  strain-displacement  relations  (2.7),  equations  (2.10), 
with  C  >  0,  represent  a  system  of  equations  to  determine  Che  unknowns 
u,  5,  0,  These  equations  must  be  solved  subject  to  certain  initial 
and  boundary  conditions.  Here,  we  define  the  initial  conditions  at  each 
point  of  Che  region  P  by  specifying 


u-u(0“)  ,  6-6(0“)  ,  0  -  0^(0“)  ,  0  -  *  (0“)  ,  (3.4a,b,c,d) 

^  fi  11 


u  -  V  (0“)  ,  6  -  w  (0“)  ,  at  t  -  t 

~  -.<0'  '  ’  r 


(3.Ae,f) 


where  u  ,  v  ,  6  ,  w  ,  0  ,  ^  are  specified  functions  of  0“  only.  Mixed 
~o  ~o  ~o  o  o  ^ 

boundary  conditions  at  each  point  s  of  the  boundary  6P  may  be  defined 
by  specifying  either: 


n  -  n(s,t)  or  u  -  u(s,t)  , 


(3.5a,b) 


m  -  m(s,t)  or  6  -  6(s,t)  , 


(3.5c,d) 


k  -  k(s,t)  or  0  -  0(s,t)  , 


(3.5e,f) 


Mixed-mixed  boundary  conditions  may  be  specified  but  are  not 
considered  explicitly. 

*The  temporary  use  of  the  symbol  s  for  a  point  on  aP  should  not  be 
confused  with  Che  use  of  the  same  symbol  elsewhere  for  the  external 
entropy  supply. 


-  kj(8,t)  or 


(|»(8,t)  ,  for  t  E[t^,  «) 


(3.5g.h) 


More  general  boundary  condltlona  than  (3.5),  which  Include  mechanical 
contact  with  a  linear  elaatlc  media  or  thermal  radiation,  ray  be  written 
at  each  point  8  of  QP  in  the  generalized  form* 


n  +  B(8)  u  +  C(a,t)  •0  ,  m  +  B  (a)  6  +  C  (8,t)  *  0  ,  (3.6a,b) 


k  -  B(8,t)(9  -  0)  +  C(8,t)  -  0  ,  -  B^(8,t)  (|)  +  C^(8,t)  -  0  , 


(3.6c,d) 


where  B  and  B  are  aaeumed  to  be  contlnuoua  aymmetrlc  three-dlmenalonal 

tenaor  functions  of  a  and  are  Independent  of  time;  B  and  B^^  are 

continuous  scalar  functions  of  (8,t);  C  and  C,  are  continuous  vector 

~  ~1 

functions  of  (a,t);  and  C  and  are  continuous  scalar  functions  of 
(s,t).  For  later  convenience,  we  define  the  scalars 


J,  -  u  •  B  u 


7  6  •  B,  6  . 


(3.7a,b) 


A  boundary  condition  of  the  type  (3.5a)  can  be  obtained  trivially  from 
(3.6a)  by  setting  B  ■  0,  and  C  ■  -  n(s,t),  and  a  boundary  condition  of 
the  type  (3.5b)  can  be  obtained  from  (3.6b)  by  setting  the  tensor  B  equal 
to  a  scalar  b  times  the  Identity  tensor  i(B  «  b  I)  and  C  -  b  u(8.t), 
and  then  taking  the  limit  as  b  approaches  Infinity. 

Recall  from  [4]  that  the  assigned  fields  f,  X,  s,  and  s^  Include 
contributions  from  both  the  effects  of  three-dimensional  body  force  and 
external  entropy  supply  as  well  as  from  the  effects  of  surface  tractions 
and  entropy  flux  on  the  major  surfaces  of  the  shell.  In  view  of  the 
specification  9^  *  ±  h/2  defining  the  major  surfaces,  we  may  write  these 
assigned  fields  In  the  forms 


A  condition  of  the  type  (3.6c)  has  previously  been  considered  for  the 
three-dimensional  theory  [8,  Sec.  5.7]. 


»  -*^  • 
l'> 


Xf 

XA 


Xs 


Xs 


1 


Xf  +  A^'^^(b^  t'^  +  b“  t”)  , 

(3.8a) 

.  .l/2/h\,  +  +  —  — 

XA  +  A  (j)  (b  t  -  b  t  )  , 

(3.8b) 

Xs  -  A^'^^(b'^  k"*”  +  b"  k")  , 

(3.8c) 

Xi^  -  A^^^(|)(b'^  k^  -  b"  k“)  , 

(3.8d) 

where  f.  A,  s,  s^  are  considered  to  be  specified  assigned  fields 
associated  with  the  three-dimensional  body  force  and  external  supply  of 
entropy,  b^  and  b~  are  positive  scalar  functions  of  the  variables  (2.17) 
and  are  Independent  of  time,  t^  la  the  traction  vector  and  Is  the 

entropy  flux  on  the  major  surface  5P'*',  and  t~  is  the  traction  vector  and 
k”  Is  the  entropy  flux  on  the  major  surface  5P~.  To  allow  mechanical 
contact  of  the  major  surfaces  with  an  elastic  media  and  to  allow  thermal 

radiation,  we  assume  that  at  each  point  of*  bP'*’ 

t  +  B  u  +  C  -  0  ,  k  -  B(e  -®)+C  =  0  (3.9a,b) 

and  at  each  point  of  bP' 

t"  +  b”  u~  +  c”  -  0  ,  k~  -  B~(9~  -  0)  +  C~  -  0  ,  (3.ina,b) 

where  b^  and  B  are  continuous,  symmetric,  three-dimensional  second 
order  tensor  functions  of  the  variables  V  In  (2.17)  and  are  Independent 
of  time;  and  C  are  continuous  vector  functions  of  the  variables 
(V,t);  and  B"*",  B  ,  C^,  C~  are  continuous  scalar  functions  of  the 
variables  (V,t).  For  later  convenience,  we  define  the  scalars 


Note  that  the  range  of  the  convected  coordinates  9“  on  the  major 
surfaces  bP"*”  and  bP“  Is  the  same  as  that  on  the  reference  surface  P  of 
the  Cosserat  surface. 
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■  AV-'A'.--  ’ 


(3.11a,b) 


J3  ”  Y  “  •  S  Ji 


J,  “  ^  u 
A  2  ~ 


Depending  on  the  nature  of  the  problem  to  be  considered,  the  quantities 
t^,  t  ,  k^,  k  ,  u^,  u  ,  6^,  9  are  either  specified  or  to  be  determined 
by  the  field  equations.  For  example,  if  we  were  to  consider  a  slightly 
more  general  boundary  condition  of  the  type  (3.9a),  which  specified  the 
tangential  components  of  t^  and  the  normal  component  of  u^,  then  the 
normal  component  of  t  and  the  tangential  components  of  u  would  be 
determined  by  the  field  equations.  This  is  similar  to  the  use  of  such 
quantities  in  the  theory  of  laminated  composite  plates  [9]. 


A.  A  Uniqueness  Theorem 

We  now  state  the  following  uniqueness  theorem:  Let  u,  6,  9,  4)  be 
displacements  and  temperature  fields  that  satisfy  the  above-mentioned 
linear  field  equations,  constitutive  equations,  and  statement  of  the 
second  law  of  thermodynamics  on  P  x  [t^,  *),  and  satisfy  the  initial 
conditions  on  P  at  t  ■  tg,  the  boundary  conditions  on  ftP  x  [tg,  “), 
conditions  of  the  type  (3.9)  on  5P^  x  [tg,  ®)  and  (3.10)  on  ftP”  x  [tg,  ®) , 
for  prescribed  values  of  the  assigned  force  f ,  director  couple  I,  and 
external  supplies  of  entropy  3  and  s^.  Then,  provided  the  specific 
kinetic  energy  K,  defined  by 


K  ■  K(v,w)  (v«v+2y  v»w+y  w»w)  , 


(4.1) 


is  positive  definite,  the  specific  heats  and  in  (2.15a)  are 
positive  scalars,  the  portion  of  the  Helmholtz  free  energy  in  (2.15b) 
is  positive  seml-def Inlte,  and  the  scalars  J^,  J2»  J3>  B~ 

in  (3.6),  (3.7),  (3.9)-(3.11)  are  positive  semi-definite,  there  exists 
at  most  one  set  of  functions  u,  6,  9,  *  that  satisfies  the  strain- 
displacement  relations  (2.8),  the  field  equations  (2.10)  (with  ^  =  0  in 
2.10d)  and  (2.12),  the  constitutive  equations  (2.15)  and  (2.19),  the 
restriction  (2.20a),  initial  conditions  (3. A),  boundary  conditions 
(3.6),  and  conditions  (3.9)  and  (3.10),  are  of  class  on  OP  x  [tg,  =») , 

and  are  of  class  C  on  P  x  [tg,  <») .  For  convenience,  the  restrictions 


-  W?.'  •.« 


'■< j'*  "S  Vis'! 


stated  above  may  be  written  In  the  mathematical  forms 


-  (yb  >0  ,  P3  >  0  ,  >  0  ,  >  0  , 


(A.2a,b,c,d) 


J,  >  0  .  >  0  t  J,  >  0  ;  J,  >  0  , 

1-  2-  3-  4- 


(4.2e,f,g,h) 


B>0  ,  Bj>0  ,  B>0  ,  B>0  . 


(4.21, j.k,^) 


Apart  from  the  discussion  of  the  generalized  boundary  conditions, 
our  method  of  proof  Is  nearly  Identical  to  that  used  In  Reference  [6] . 
Specifically,  we  assume  the  existence  of  two  different  solutions  of  the 
abovu-stated  Initial,  mixed  boundary-value  problem,  form  the  difference 
solution,  and  use  a  consequence  of  the  field  equations  to  prove  that  the 
difference  solution  Is  the  null  solution.  Let  us  denote  the  two 
solutions  by  the  sets  of  quantities 

-f  —  +  - 

U  -  {u,  6,  0,  ♦,  n,  k,  m,  Ti,  q, ,  ,  p,  p, ,  t  ,  t  ,  k  ,  k  }  (4.3a) 

^  ^  ^  ^  ^  J,  X  ^  ^JL  ^  ^ 

U'  -  {u'.  6'.  0'.  £'.  k',  o’,  h’,  Ti{.  CJ,  e',  £’■*■,  t’”,  k’'*',  k’”} 

(4.3b) 


and  form  the  difference  solution 


U  -  U  -  U’  . 


(4.4) 


It  follows  that  the  difference  solution  satisfies  the  following:  the 
field  equations 


\(v  +  y^  w)  -  A^^^(b'’'  t'*'  +  b”  t")  +  (A^^^  n“) 


(4.5a) 


The  restriction  (4.2b)  Is  consistent  with  the  condition  (2.20b). 


\(y^  V  +  w)  -  -  b"  t')  -  k  +  (A^'*^  M®) 


1/2  'o. 


XtI  -  -A^^^Cb"^  k'*'  +  b"  k")  -  (A^^2  p“) 

»a 


1  -  A^/2(|](b+  k+  -  b“  k")  -  (A^'^  ;®)^^  , 


1/2  -a. 


(4.5b) 


(4.5c) 


(4.5d) 


on  P  X  [t^,  »),  the  restriction 


Po  ♦?  +  £*£  +  £i  *  Jli  < 


on  P  X  [t  ,  «),  the  initial  conditions 


(4.6) 


u*0  ,  6*0  ,  9«0  ,  (>“0  ,  u“0  ,  6*0  (4.7a-f) 


on  P  at  t  ■  t  ,  the  boundary  conditions 


n  +  Bu*0  ,  m  +  B,  6*0  , 


k-B0-O  ,  kj-Bj<|»»0  , 


on  dP  X  [t  ,  •),  the  conditions 


t'^  +  b'*’  u'*’  -  0  ,  k"^  -  b"*"  e'*'  -  0  , 


on  91P^  X  [t  ,  •),  and  the  conditions 


t+Bu-O,k-B0-O, 


on  6P'  X  [tg,  «). 


(4.8a,b) 


(4 .8c,d) 


(4.9a,b) 


(4.10a,b) 


Multiplying  (4.5c)  by  0,  (4.5d)  by  adding  the  results  together, 
integrating  over  the  region  P,  using  the  divergence  theorem  and  the 
conditions  (4.8c,d),  (4.9b),  and  (4.10b),  we  obtain 


*«■  «-«  •••  «•  ••  ••  . 

v '.'.V  ‘v  ‘v  V  ■-•  ■ 


^  •".'A . 


.0*- “  •  v.v  • • 


(4.11) 


/p  PqO  n  +  ♦  Ti^)da  -  /p(p^  ♦  +  £  •  g  +  •  ^{>da 

/>  2  ^2 

-  /p[t>^  b'^(9'^)  +  b"  (9")  ]da 

-  /^p(B  +  Bj  i^)d8  , 

where  do  Is  the  area  element  on  P  and  ds  Is  the  arc  length  on  dP.  Taking 

^  A 

the  Inner  product  of  (4.5a)  with  v,  (4.5b)  with  w,  adding  the  results 
together.  Integrating  over  the  region  P,  using  the  divergence  theorem  and 
the  conditions  (4.8a,b),  (4.9a),  and  (4.10a),  we  deduce  the  expression 

•  • 

E  -  /p  p^(e  ^|  +  i  ipdo  (4.12) 


where 

E  -  /p  [po(4>i  +  K)  +  y  P3  0^  +  Y  +  b"^  J3  +  b’  J^]do 

+  /^p(J^  +  J2)ds  ,  (4.13a) 

A  AAA  AAA  A  A 

’•'1  “  '•'l^®ap’  ^1*  ^^la^  ’  ^  »  (4.13W) 

J2  -  ^2^^^  *  '^3  "  '^3^'"  >  »  *^4  •  (4.13e,f,g) 

and  where  the  functions  4»j,  K,  -  J^,  are  defined  by  (2.15b),  (4.1), 
(3.7),  and  (3.11),  respectively.  Now  from  (4.11)-(4.13)  and  the 
restrictions  (4.2),  we  realize  that  E  >  0  and  E  <  0.  Using  the  Initial 
conditions  (4.7),  we  obtain  the  result  that  E  -  0  for  all  time  t  cft^,  *) 
and  therefore 

A  A  A  A 

v*0  ,  w*0  ,  9»0  ,  ^“0  .  (4.14a,b,c  ,d) 

Integrating  (4.14a,b)  and  using  the  initial  condition  (4.7),  we  have 

A  A 

u  -  0  ,  6-0  ,  (4.15a,b) 
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Finally,  substituting  (4.14c,d)  and  (4.15a,b)  Into  the  conditions  (4.9), 

(4.10) ,  we  conclude  that 

t^-0  ,  t  -0  ,  k^-0  ,  k  -0  (4.16a,b,c,d) 

which  completes  the  proof. 

To  prove  uniqueness  for  thermoelastic  statics,  we  need  slightly 
stronger  conditions  than  (2.20a)  and  (4. 2d).  Specifically,  we  retain 
the  restrictions  (4.2e-jl)  and  assume  that  (|>|^  Is  positive  definite  and 
that  the  expression  on  the  left-hand  side  of  (2.20a)  Is  negative 
definite  so  that 

4,^  >  0  ,  Pq  ♦  +  £  •  S  +  <  0  (4.17a,b) 

where  In  (4.17a)  vanishes  only  when  the  mechanical  fields  Yi> 

vanish  and  (4.17b)  vanishes  only  when  the  thermal  fields  9 
^  -  vanish.  Furthermore,  we  require  the  temperature  9  to  be  specified 

t 

on  at  least  one  point  on  the  boundary  of  the  shell.  This  can  be  done  by 
specifying  9  on  8P,  9"*"  on  bP'*’,  or  9”  on  SP". 

Now,  for  thermoelastlc  statics  the  thermal  equations  (4.5c,d)  are 
uncoupled  from  the  mechanical  equations  (4.5a,b)  and  the  expression 

(4.11)  can  be  derived  with  the  left  hand  side  vanishing.  It  follows 

^  * 
that 

Po  ♦  ^1  El  • 

from  which  we  conclude  that 

A  A 

g  "  0  ,  ♦  -  0  . 


*Recall  from  (2.15f)  that 


(4.19a,b) 


Integrating  (4.19a)  and  using  the  specification  of  6  at  a  point  of  the 
boundary  of  the  shell,  we  have 


e  -  0  .  (4.20) 

Substituting  (4.19b)  and  (4.20)  Into  the  conditions  (4.9b)  and  (4.10b), 
we  deduce  that 

it'’"  -  0  ,  k"  -  0  ,  (4.21a,b) 

which  completes  the  proof  for  the  thermal  fields.  In  view  of  the 
results  (4.19b)  and  (4.20),  we  may  take  the  Inner  product  of  the 

A 

equilibrium  form  of  equations  (4.5a)  and  (4.5b)  with  u  and  5,  respec¬ 
tively,  and  derive  the  expression 


/p  (Pq  \  J3  +  b“  J^)d<j  +  /^p  (J3  +  J2)ds  -  0  .  (4.22) 

A 

It  follows  that  (|f^  a  0  so  that 

>  ■  "  '  n  ■  "  •  'la  ■  O  •  (4.23., b,c) 

Hence,  the  displacements  are  unique  to  within  a  linear  superposed  rigid 
body  displacement.  If  this  arbitrariness  Is  removed,  then  the  displace¬ 
ments  will  be  unique 

A  A 

u  ■  0  ,  6*0  (4.24a,b) 

and  from  the  conditions  (4.9a)  and  (4.10a),  we  can  conclude  that 
«_ 

t  ■  0  ,  t  a  0  ,  (4.25a,b) 

which  completes  the  proof  for  the  mechanical  fields. 
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Appendix  B 

A  NONLINEAR  CONSTRAINED  THEORY  OF  SHELLS  THAT 
INCLUDES  TANGENTIAL  SHEAR  DEFORMATION 


(Submitted  for  publication  to  ASME  Journal 
of  Applied  Mechanics) 


Introduction 


Within  the  context  of  theories  of  deformable  bodies,  a  constrained 
material  Is  a  material  that  can  experience  only  a  restricted  class  of 
motions.  Various  constrained  theories  of  shells  have  been  developed 
mainly  because  the  system  of  equations  characterizing  a  constrained 
theory  Is  simpler  than  that  characterizing  the  general  theory.  For 
shells.  It  Is  common  to  develop  constrained  theories  that  exclude  one  or 
both  of  the  following  two  types  of  deformations:  (a)  normal  extension 
and  (b)  tangential  shear  deformation.  The  terminology  normal  extension 
and  tangential  shear  deformation  is  used  Instead  of  the  usual  termi¬ 
nology  "transverse  normal  strain,"  and  "transverse  shear  deformation" 
because  It  Is  more  descriptive  when  considering  nonlinear  deformation  of 
a  shell  as  opposed  to  linear  deformation  of  a  plate. 

To  discuss  constrained  theories  of  shells  of  the  type  considered 
here.  It  Is  particularly  convenient  to  model  the  shell  as  a  Cosserat 
surface  [1].  Although  various  restricted  or  constrained  theories  of 
shells  [I,  Sec.  10;  2,3,4]  have  been  developed,  we  choose  to  focus 
attention  on  three  constrained  theories.  For  constrained  theory  I,  we 
exclude  both  deformations  (a)  and  (b)  so  that  a  material  fiber  that  Is 
Initially  normal  to  the  undeformed  reference  surface  of  the  shell 
remains  normal  to  the  deformed  reference  surface  and  has  constant 
length.  A  nonlinear  theory  of  this  type  has  been  developed  as  a 
restricted  theory  [1,  Sec.  10],  and  it  has  been  observed  that  It  reduces 
to  the  classical  linear  Klrchhoff-Love  plate  theory.  For  constrained 
theory  II,  we  exclude  only  deformation  (b)  so  that  a  material  fiber  that 
Is  Initially  normal  to  the  undeformed  reference  surface  of  the  shell 
remains  normal  to  the  deformed  reference  surface,  but  Is  allowed  to 
extend  or  contract  In  length.  A  nonlinear  theory  of  this  type  has  also 
been  developed  [2].  For  constrained  theory  III,  we  exclude  only 
deformation  (a)  so  that  a  material  fiber  that  Is  Initially  normal  to  the 
undeformed  reference  surface  of  the  shell  Is  allowed  to  deform  away  from 
the  normal  to  the  deformed  reference  surface,  but  the  component  of  the 
material  fiber  normal  to  the  deformed  reference  surface  remains  constant. 
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Although  a  theory  of  this  type  can  be  reduced  to  a  linearized  theory  of 
a  plate  that  Includes  "transverse"  shear  deformation  [5],  a  nonlinear 
version  of  this  theory  does  not  appear  to  have  been  previously  developed. 


Our  purpose  here  Is  to  develop  a  nonlinear  version  of  constrained 
theory  III.  Because  the  constraint  associated  with  this  theory  Is 
mechanical  In  nature,  we  will  confine  attention  to  the  purely  mechanical 
theory.*  In  the  following  sections,  we  record  the  basic  equations  and 
discuss  the  constraint  and  constraint  response  associated  with  theory 
III.  Next,  the  boundary  conditions  are  discussed,  with  a  few  comments 
on  the  linear  theory,  and  finally  the  Initial  conditions  are  stated. 

Basic  Equations 

In  this  section,  we  use  the  notation  In  [2]  and  briefly  record  the 
basic  equations  appropriate  for  a  constrained  version  of  the  theory  of  a 
Coaserat  surface.  For  a  complete  discussion  of  this  theory,  we  refer 
the  reader  to  [1,6].  Let  the  material  points  of  the  Cosserat  surface  C 
be  Identified  by  means  of  a  system  of  convected  coordinates  0®(a  ■  1,2) 
and  let  the  two-dimensional  region  occupied  by  the  material  surface  in 
the  present  configuration  at  time  t  be  denoted  by  c.  Further,  with 
reference  to  the  present  configuration,  let  the  vector  valued  function 
r  define  the  position  of  a  material  point  of  C  and  at  each  such  point 
define  the  vector  valued  function  d,  called  the  director.  Then,  a 
motion  of  the  Cosserat  surface  Is  defined  by 

r  -  r(9®,t)  ,  d  -  d(9®,t)  ,  [a^  £2  ^  (la.b,c) 

where  the  tangent  vectors  a  and  the  unit  normal  vector  a^  are  defined 
by 


*A  thermomechanical  theory  of  the  type  [6]  could  be  developed  by 
appropriately  modifying  the  constitutive  equations. 


a,  •  a,  -  1 


r 


i 


a  -  -  ,  a  •  a-  -  0  , 

~o  Qgtt  ’  ~o  ~3 


.  a"'^  -  fa^  £2  £3]  >  0 


(2a,b,c,d) 


The  velocity  v  and  director  velocity  w  are  defined  by 


V  ■  r  ,  w  ■  d  , 


(3a, b) 


where  a  superposed  dot  denotes  material  differentiation  holding  9°^ 
fixed . 

In  the  reference  configuration,  the  various  kinematic  quantities 
may  be  denoted  by* 


r-R  ,  d*D*DA,  ,  a.  *A, 


1/2  ,1/2  ..  . 

a  -  A  ,  (4a,b,c,d) 


1  /  ^  /r 

where  R,  D,  A^^,  A  depend  on  tlM  coordinates  9“  only.  The 
specification  (4b)  la  made  without  loss  In  generality.  For  later 
convenience,  we  recall  [1]  the  kinematic  definitions 


a  -  •  a  'a.  ,  d,  •  a,  •  d  ,  -  a.  •  d 

aP  ~o  ~p  1  ~1  ~  la  ~1  ~,a 


(5a,b,c) 


®ap  “  2  ^*ap  -  ^p>  »  Yi  -  (5d,e,f) 


b  .  -  a  „  •  a,  , 
aP  ~a,P  ~3 


where  A^^,  D^,  are  the  reference  values  of  a^^,  d^,  b^^, 

respectively;  e  y.,  k  are  strains  measured  relative  to  the 
aP  1  la 

reference  configuration;  and  a  comma  denotes  partial  differentiation 
with  respect  to  9®. 


Throughout  this  text,  Greek  Indices  have  the  range  (1,2)  and  Latin 
Indices  have  the  range  (1,2,3). 


Let  P,  bounded  by  the  closed  curve  bP,  denote  the  region  occupied 
by  an  arbitrary  material  portion  of  the  surface  c  In  the  present 
configuration  and  let  v  be  the  unit  outward  normal  to  9P.  Using  the 
notation  of  [2]  and  referring  to  the  present  configuration,  we  define 
the  following  quantities:  the  mass  density  (mass  per  unit  area  of  P) 
p  ■  p(9*,t);  the  contact  force  n  -  nCQ^'.t;  v)  and  the  contact  director 
couple  m  -  m(9*,t;  v),  each  per  unit  length  of  the  curve  dP;  the 
specific  (per  unit  mass  of  P)  assigned  force  f  «  f(9**,t)  and  specific 
assigned  director  couple  i,  «  1(9®, t);  the  Intrinsic  director  couple 
k  ■  k(9®,t)  per  unit  area  of  P;  and  the  Inertia  coefficients  y^  *  y^(9“) 
and  y^  ■  y^(9®),  which  are  Independent  of  time.  With  suitable  continuity 
assumptions.  It  can  be  shown  [1]  that 


n»N  V  ,  m*Mv  ,v“a  ‘v  , 
o.  ~  a  ~  ~  a  s  ^ 

.a  .  .A 


(6a,b,c) 


where  N  and  M  are  Independent  of  v.  Referring  all  tensor  quantities 

.  ,  ** 
to  base  vectors  a  ,  we  may  write 


„a  „ia  ,  ,  1  „la 

N*N  a.  ,k*ka,  ,M«M  a.  , 


(7a,b,c) 


where  the  usual  summation  convention  Is  used. 

Coupled  thermomechanical  constraints  for  a  three-dimensional 
continuum  have  been  previously  discussed  within  the  context  of  more 
classical  thermodynamics  [7]  as  well  as  within  the  context  of  recent 
developments  [8].  Further,  a  rather  general  discussion  of  purely 
mechanical  constraints  for  the  theory  of  a  Cosserat  surface  Is  contained 


*Thls  notation  differs  from  that  used  In  [1,6].  In  particular,  we  note 
that  the  quantities  a,  Jj,  jj  defined  here,  correspond, 
respectively  to  fi,  jg,  JJ  In  [1,6]. 

The  notations  used  here  for  and  M®  are  consistent  with 
those  of  [2],  but  differ  from  the  notation  and  Ij“  » 

used  In  [1,6]. 


In  [9].  In  these  works,  a  constrained  material  Is  characterized  by  a 
set  of  constraint  equations  that  restrict  various  klnematlcal  quantities 
(deformation  and  temperature  fields)  as  well  as  associated  constraint 
responses  that  Introduce  a  certain  arbitrariness  In  the  kinetic 
quantities. 

For  constrained  theory  III,  we  need  Introduce  only  a  single 
constraint,  which  excludes  normal  extension  and  Is  characterized  by 


d»a-»D«A-»n 


(Y3  *  0  .  Y3  ="  0) 


Following  previous  works  [6-8],  we  assume  that  the  kinetic  quantities 
^oc6  flt  iot 

N  k  ,  M  are  completely  determined  by  constitutive  equations  and 
3 

that  k  separates  Into  two  additive  parts:  one  part,  denoted  by  k  ,  Is 

-3 

determined  by  a  constitutive  equation;  and  the  other  part,  k  ,  Is  an 
arbitrary  function  of  (0®,t),  Independent  of  strain  rates,  which  Is 
further  assumed  to  be  workless.  Thus,  we  assume  that 


3  *3  -3  -3  • 

k  -  k^  +  k^  ,  k  Y,  -  0  . 


(9a, b) 


In  view  of  the  constraint  (8),  the  director  d  may  be  expressed  In  the 


d  ■  d  3  +  D  3m 


It  Is  now  convenient  to  record  [1,2]  the  local  forms  of  the  basic 
equations  of  motion  as: 


^  .1/2 
\  ■  p  a  ■  A 


(11a) 


p  f“  +  N“P|p  -  b®  -  0  ,  p  f”'  +  +  N®P  '’ap  -  0  . 


(11b. c) 


*The  quantity  will  be  defined  presently. 
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p  r  -  k-  +  bp  -  0 


(lid) 


(lie) 


where 


n“P  -  nP“  -  n“P  -  X?a  -  d“(p  e  +  m“”|  -  bP) 


P  ^  uPl 


(12a) 


.3a 


»,®P|  wPP  J®/”  13  ,  ,3P|  .  ap  ,  ■(  .  ,3  aa  ,a  .3a 

d  ( p  -t  +  M  |.  -  M  b„ )  —  d(pJl  +  M  I-+M  ^_o)  ■*"  ^  ~  ^ 


aP' 


(12b) 


*1  “1  i 

and  where  f  and  I  are  defined  using  the  vectors  a  which  are  the 

reciprocals  of  a  ,  such  that 

^1  1  •  1  •  ^1  1  1  •  2  • 
f  -a  ‘(f-v-y  w)  ,  Jl  -  a  •  (I  -  y  v-y  w)  .  (13a, b) 

In  (11a)  Pq  Is  the  reference  value  of  the  mass  density  p.  Also,  in  (11) 
and  (12)  a  bar  denotes  covariant  differentiation  with  respect  to  the 
metric  a^p.  Equation  (11a)  represents  conservation  of  mass,  equations 
(11b, c)  represent  the  balance  of  linear  momentum,  equations  (lld,e) 
represent  the  balance  of  director  momentum,  and  equations  (12a, b) 
represent  the  balance  of  angular  momentum. 

Once  appropriate  constitutive  equations  are  specified  for  the 
kinetic  quantities 

.  k*.  M**)  ,  (14) 


the  six  equations  (lla-d)  may  be  used  to  determine  the  six  unknowns 

{p.  r,  d®}  .  (15) 

*3 

For  an  arbitrary  value  of  k  ,  equation  (lie)  may  be  satisfied  by  an 

-3 

appropriate  specification  of  k  .  Now,  equations  (llb-d)  must  be  solved 


subject  to  certain  boundary  and  initial  conditions,  which  are  discussed 
in  Che  next  section. 


Boundary  and  Initial  Conditions 


A  rather  detailed  discussion  of  boundary  conditions  for  constrained 
theories  I  and  IT  is  Included  in  [1,  Sec.  15]  and  [2],  respectively. 
Basically,  we  follow  the  discussion  in  [1]  and  recall  [2]  that  the 
boundary  integral  appearing  in  the  energy  equation  cakes  the  form 

/  (n  •  V  +  m  •  w)ds  »  /  (n^  v.  +  w. )ds  ,  (16) 

bP  ~  ~  ~  ~  ap 

where  ds  is  the  elemental  arclength  of  dP  and  where 


1  1  „la  i  t  ^la 

n  “a  ‘n^N  v  ,  m“a  v 

~  ~  a  ’  .w  a 


V,  *  a,  •  V  ,  w.  ■  w  •  a^ 
i  ~1  ~  1  ~  ~1 


(17a, b,c) 


(17d,e) 


Recalling  [1],  that  an  arbitrary  function  F  »  F(9“,t)  may  be  expressed 
in  terms  of  its  normal  and  tangential  derivatives  dF/av  and  aF/as, 
respectively,  we  may  integrate  (16)  by  parts,  assume  continuity  on  aP, 
and  use  (10)  to  rewrite  (16)  in  the  form 


,  ov  ov 

/  [P  (v, )  +  m'^(3  -  V  d^  -  V  D  -r — ) 

^  i  ^  a  a  av  o  av 

1  o 

''p 

where  P^  are  given  by 

P"^  -  [n'^  +  m^(d^  -  D  bp  +  m^  d^  bp  +  [m^  d*’  ,  (19a) 


P^  -  (n^  +  m^  d*^  b^)  +  [(m^  P  -  m^  * 


(19b) 


and  where  the  unit  tangent  vector  X  to  aP  and  Che  Chrlstoffel  symbol 
rf-  are  defined  by 

\p 


It  Is  clear  from  (18)  that  at  each  point  of  3F,  we  must  specify  either 
the  kinetic  variables  (P^,  m^)  or  their  associated  kinematic  variables 
that  are  written  In  parentheses  In  (18).  Thus,  for  the  nonlinear 
constrained  theory  III,  we  must  specify  six  boundary  conditions  at  each 
point  of  bP.  However,  for  a  theory  that  Is  linearized  about  the 
reference  configuration,  the  quantities  P^,  m^,  d^,  v^  are  small  and  of 
order  e(o  <  e  «  1)  so  that  the  term  In  (18)  associated  with  m^  Is 
negligible  compared  to  the  other  terms.  Thus,  for  the  linearized  theory 
the  number  of  boundary  conditions  reduces  to  five. 

To  see  that  the  terms  In  (18)  associated  with  bending  of  a  plate 
are  consistent  with  those  discussed  In  [5],  we  note  that  for  the  linear 
theory  of  a  plate 


o  .,3  3  ,  a  ,  <j  . 

P  -  n  ,  P  -  n  +  ^  (m  D)  , 

bv*  bv- 

w  -  (d  -  D  V,  )a‘’  -  (5  -  D  T-^IX  +  fd  -  D 
~  'a  3,a  ~  '■  s  bs  v  bv 


(21a, b) 
(21c) 


d  “X’d  ,  d  “v^d  . 


(21d,e) 


where  we  have  Introduced  the  temporary  notation  dg  and  d^,  respectively, 
for  the  tangential  and  normal  components  of  d  relative  to  the  curve  bP. 
Thus,  for  a  plate,  (18)  may  be  rewritten  as 


.  »  bv-  bv, 

/  [n^  V  +  -2-  (m'^  X  D  v,)  +  m*^  (3  -  D  X  -  D  v  ^l]ds 
gp"-  lbs'  a  3'  '■o  obs  obv^-* 


(22) 


Integrating  the  second  terra  In  (22)  and  using  continuity,  we  may  rewrite 
(22)  In  the  form 


/  [(n  •  X)(v  •  X)  +  (n  •  v)(v  •  v)  +  n  v_ 
bP 

bv-  bv- 

+  <5  ‘  ^)(^s  ■  °  bi^^  ° 
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Next,  recall  (1]  that  the  Cosserat  theory  may  be  brought  Into  a  one-to- 
one  correspondence  with  the  three-dimensional  theory  If  the  position 

A 

vector  £  locating  points  in  the  shell  admits  the  representation 

£*  -  r  +  d  (24) 

where  9  Is  the  convected  coordinate  through  the  thickness  of  the  shell. 

Using  the  notation  u^  and  u^  defined  In  [5],  It ^follows  from  (21c)  and 

(24)  that  w  •  \  and  w  •  v  correspond  to  u'  and  u' .  respectively,  and 
~  ~  ~  ~  s  n 

that  the  last  three  terms  In  (23)  correspond  to  a  dynamical  version  of 
the  boundary  Integral  In  [5],  where  we  note  that  D  -  1  characterizes  a 
shell  of  constant  thickness. 

Finally,  for  a  dynamical  problem  It  Is  necessary  to  specify  Initial 
conditions  of  the  form 

£  *  £(®“)  •  X  *  »  S  "  »  for  t  -  0 

(25a,b,cd) 

where  r,  v,  d,  w  are  specified  functions  of  0®  only  and  d,  w  are 
consistent  with  the  constrained  form  (10)  of  the  director. 
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Appendix  C 


ON  THE  DETERMINATION  OF  CERTAIN  CONSTITUTIVE 
COEFFICIENTS  FOR  THERMOELASTIC  SHELLS 


(Submitted  for  publication  to  the  International 
Journal  of  Solids  and  Structures) 


ABSTRACT 


In  this  paper  we  model  a  shell  composed  of  a  linear  elastic, 
homogeneous,  isotropic  material  as  a  Cosserat  surface.  Specific 
attention  is  focused  on  the  determination  of  certain  thermal  and 
mechanical  constitutive  coefficients,  which  were  previously  determined 
by  integrating  three-dimensional  constitutive  equations.  Here,  we 
determine  these  coefficients  by  comparing  Cosserat  solutions  with  exact 
three-dimensional  solutions.  This  comparison  suggests  values  for  one  of 
the  thermal  coefficients  and  two  of  the  mechanical  coefficients  that  are 
different  from  those  previously  proposed  [1,2]. 


C-2 


1 


Introduction 


Classical  developments  of  shell  theory  usually  start  with  a 
complete  three-dimensional  description  of  a  shell-llke  body  Including 
constitutive  equations.  Shell  equations  that  are  functions  of  two 
space- variables  and  time  are  then  developed  by  Introducing  approxima¬ 
tions  and  either  Integrating  the  equations  of  motion  through  the 
thickness  of  the  shell  or  by  using  an  Integral  form  of  a  variational 
principle.  More  recently.  It  has  become  common  to  develop  shell 
equations  by  a  direct  approach  In  which  the  shell  Is  modeled  as  a 
Cosserat  surface  [1].  This  latter  approach  has  the  distinct  advantage 
over  the  classical  approach  that  equations  can  be  developed  for  shells 
with  arbitrary  constitutive  properties.  More  specifically,  the 
discussion  of  constitutive  equations  within  the  context  of  the  Cosserat 
theory  Is  very  similar  to  that  In  the  three-dimensional  theory. 
Constitutive  coefficients  are  determined  by  comparing  the  predictions  of 
shell  theory  with  experimental  data  or  exact  solutions  of  the  three- 
dimensional  equations  that  have  already  been  shown  to  accurately  predict 
experimental  data. 

Until  recently,  the  Cosserat  theory  admitted  any  finite  number  of 
directors  to  provide  Information  about  the  variation  of  mechanical 
variables  through  the  thickness  of  the  shell  but  only  admitted  a  single 
temperature  field  to  model  the  average  temperature  In  the  shell.  Recent 
advances  In  thermodynamics  have  allowed  the  Cosserat  theory  to  be 
generalized  to  admit  any  finite  number  of  temperature  fields  that 
provide  Information  about  the  variation  of  temperature  through  the 
thickness  of  the  shell  [2].  An  Important  special  case  of  the  general 
theory  Is  one  that  admits  a  single  director  and  two  temperature  fields. 

In  [2]  specific  constitutive  equations  are  developed  for  a  plate 
composed  of  a  linear  elastic,  homogeneous.  Isotropic  material.  Values 
for  most  of  the  mechanical  constitutive  coefficients  and  for  some  of  the 
thermal  coefficients  have  been  determined  by  comparing  predictions  of 
the  Cosserat  theory  with  exact  solutions  of  the  three-dimensional 
equations  [1].  Values  for  the  mechanical  constitutive  coefficients 


associated  with  tangential  shear  deformation  (more  commonly  referred  to 
as  transverse  shear  deformation)  and  the  new  thermal  coefficients  were 
determined  by  Integrating  three-dimensional  constitutive  equations 
[1,2].  Even  though  a  one-to-one  correspondence  between  the  three- 
dimensional  theory  and  the  Cosserat  theory  may  be  established,  these  two 
methods  of  determining  constitutive  coefficients  do  not  always  yield  the 
same  results  (see  additional  comments  In  Sections  4  and  5). 

The  objective  of  this  paper  Is  to  show  that  the  new  thermal 
coefficients  and  the  mechanical  coefficients  associated  «rlth  tangential 
shear  deformation  may  be  determined  by  direct  comparison  with  exact 
three-dimensional  solutions.  In  the  following  sections,  we  record  the 
basic  equations  valid  for  a  shell  composed  of  a  linear  elastic, 
homogeneous.  Isotropic  material.  Next,  we  briefly  recall  [1]  how  most 
of  the  mechanical  coefficients  were  determined  by  comparison  with  exact 
solutions.  Then,  we  determine  the  new  thermal  coefficients  by  comparing 
with  simple  solutions  of  the  three-dimensional  heat  conduction  equation 
for  a  rigid  plate.  One  of  the  thermal  coefficients  determined  In  this 
manner  has  a  value  different  from  that  proposed  In  [2].  Finally,  we 
determine  the  two  mechanical  coefficients  associated  with  tangential 
shear  deformation  by  comparing  with  the  exact  solutions  of  simple  shear 
of  a  plate  and  twisting  of  a  circular  cylindrical  shell.  The  values  of 
both  of  these  coefficients  determined  In  this  manner  are  different  from 
those  proposed  In  [Ij. 

2.  Basic  Equations 

Let  the  material  points  of  the  Cosserat  surface  C  be  Identified  by 
means  of  a  system  of  convected  coordinate  9*^(a  ■  1,2)  and  let  the  two- 
dimensional  region  of  space  occupied  by  the  material  surface  In  the 
present  configuration  at  time  t  be  denoted  by  c.  Further,  let  the 
vector  valued  function  x  define  the  position  of  a  material  point  of  the 
surface  C  and  at  each  such  point  define  the  vector  valued  function  d, 
called  the  director,  and  the  two  temperature  fields  6  and  4*,  each 
referred  to  the  present  configuration.  Then,  a  thermomechanical  process 
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of  the  Cosserat  surface  Is  defined  by 


r-r(9“.t)  .  d.d(9“.t)  .  £2  d]  >0  . 

9  -  9(9®, t)  ,  (9  >  0)  ,  ♦  -  ♦(9®,t)  , 

where  Che  tangent  vectors  a  and  Che  unit  normal  vector  a 

~(t  ~3 

by 

1/2 

Sa  •  ^  •  23  •  “  ■  I2i  £2 
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(2.1a,b,c) 

(2.1d,e,f) 

are  defined 

(2.2a,b,c,d) 


and  Che  condition  (2.1c)  ensures  that  Che  director  is  nowhere  tangent  to 
c.  The  velocity  x  director  velocity  j;  may  be  defined  by 


V  ■  r 


“  » 


(2.3a,b) 


where  a  superposed  doc  denotes  time  differentiation  holding  9*^  fixed. 

In  Che  reference  configuration,  we  assume  that  the  shell  has  uniform 
thickness  h  and  is  at  uniform  temperature  9^.  Then,  the  reference  values 

’jig 

of  Che  various  kinematic  quantities  may  be  denoted  by 


i,  £  » 


d  -  D  - 


~3 


(2 .4a,b,c,d) 


®  "  ®o  ’  ®  ’  (2.4e,f) 

where  R,  and  depend  on  the  coordinates  9®  only. 

Let  P,  bounded  by  the  closed  curve  9P,  denote  the  region  occupied 
by  an  arbitrary  material  portion  of  the  surface  c  in  the  present 


^Throughout  Che  text  Creek  indices  have  a  range  (1,2)  and  Latin  indices 
have  a  range  (1,2,3). 


configuration  and  let  v  be  the  unit  outward  normal  to  dP.  Using  the 

<§t 

notation  of  [3]  and  referring  to  the  present  configuration,  we  define 

the  following  quantities:  the  positive  mass  density  (mass  per  unit  area 
(Z  ct 

of  P)  p  ■  p(0  ,t);  the  contact  force  n  •  n(9  ,t;  v)  and  the  contact 
director  couple  m  -  m(9**,t;  v),  each  per  unit  length  of  the  curve  &P; 
the  specific  (per  unit  mass  of  P)  assigned  force  f  ■  f(9'’^.t)  and 
specific  assigned  director  couple  A  >  1(9*^, t);  the  Intrinsic  director 
couple  k  -  k(9*,t)  per  unit  area  of  P;  the  Inertia  coefficients  “ 
y^(9*)  and  y^  ■  y^(0*)  which  are  Independent  of  time;  the  specific 
entropies  q  ■  q(9*,t)  and  the  specific  Internal  rates  of 

production  of  entropy  l  -  5(9®, t),  -  C^(9“,t),  and  -  Ej^(9“,t);  the 

entropy  fluxes  k  ■  k(9®,t;  v)  and  k,  ■  k, (9**,t;  v);  the  specific 
external  rates  of  supply  of  entropy  s  ■  s(9®,t)  and  s^^  *  8^(9®, t);  the 
specific  Internal  energy  e  ~  e(9°,t);  and  the  specific  Helmholtz  free 
energy  4>  -  4»(9“,t)  =  e  -  9q  -  ^q^^. 

For  the  linear  theory,  It  Is  convenient  to  Introduce  the  displace¬ 
ment  vector  u  and  director  displacement  ^  relative  to  the  reference 
configuration  by  the  equations 

r*R+u,d«D+5.  (2.5a,b) 

^  ^  ^  M  M  A# 

We  now  assume  that  In  the  present  configuration,  the  displacements  u,  5 
and  the  temperatures  (9  -  9^)  and  ^  are  of  order  6(0  <  6  «  1)  and 
that  quadratic  terms  In  these  quantities  may  be  neglected  relative  to 
linear  terms.  With  this  background,  we  record  [1]  expressions  for  the 
kinematic  quantities  and  for  the  strains 

the  forms 


*Thls  notation  differs  from  that  used  In  [1,2]. 

I 

The  temporary  use  of  this  symbol  for  the  small  parameter  should  not  be 
confused  with  the  use  of  the  same  symbol  elsewhere  for  the  Internal 
energy . 


A  *  -  A  'A 
ap  “-o  ~p 


B  „  -  A  a  •  A 
aP  -<r,P  ~3 


(2.6a,b) 


‘up  ■  "‘as  •  *3a  ■  ”  ■ 


V  2  'M?  ♦  %  “3  • 


Y  *6  +  u,  ^  +  B  u.  ,  Yi  ■  6,  , 

CK  CC  5yOE  0t\  3  3 


•'aP  "  Pap  ■  ®aP  ^3  *  '^3a  “  P3a  ®a  * 


(2.6c,d,) 

(2.6e) 

(2.6f,g) 

(2.6h,i) 


PaP  “  ■“3|aP  ■  ®a|p  “x  "  \  “x|p  "  ®p  “x|a  ■*■  ®oX  ®p  “3  ‘'^ajp  ’ 

(2.6j) 


P3a 


y 


3, a 


(2.6k) 


where  for  che  linear  theory  all  censor  quantities  are  referred  to  the 
base  vectors  Aj.  In  (2.6)  and  through  the  text,  we  use  the  usual 
summation  convention  over  repeated  Indices;  a  comma  denotes  partial 
differentiation  with  respect  to  9*’^,  and  a  bar  denotes  covariant 
differentiation  with  respect  to  che  metric  A^^. 

With  suitable  continuity  assumptions,  It  can  be  shown  that*  [1,2] 


n  ■  N®  V  ■  (N^*  A, )v  ,  m  ■  M*  V  ■  (M^®  A. )v  ,  (2.7a,b) 

k-jg»v»p®v^  ,  ’  (2.7c,d) 

where  v  ■  A  •  v  are  the  components  of  che  normal  vector  v  and  where 
a  ~o  ~  ~ 

W**,  M®,  £,  and  are  Independent  of  v.  Further,  with  reference  to  the 

energy  equation,  the  specific  external  rates  of  heat  supply  r  and  r^, 

and  che  heat  flux  vectors  q  and  q,  are  defined  by 


*The  notation  In  (2.7a,b)  Is  consistent  with  that  In  [3],  but  differs 
from  chat  In  [1,2]. 
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r  ■  9  s 


a "  ®  E 


> 


ai  -  ♦  £i 


(2.8a,b,c,d) 


> 


«l>  a, 


» 


Now  the  local  forms  of  Che  basic  equations  may  be  recorded  as 


.  1/2  .1/2 
X  ■  p  a  *  P  A 
o 


(2.9a) 


f“  +  N“P|p  -  Bp  -  0  ,  f^  +  N^^lp  +  n“^  B^p  -  0  ,  (2.9b,c) 

Poi“-k%M“Plp.M^PB;.0  .  p^?-.^.M^P|p.M“%„p.O  . 

(2.9d,e) 


P  ■  P 

^O  '  ^0 


■  P  'at  ‘ 


Po  \  *  Po^®l  ^1^  ■  P 


“I 

1  'a  ’ 


(2.9f,g) 


where 


.  »«■>  .  ?  +  m“^|j  .  (Z.lOa.b) 


and  where 


“11  “11  1»2» 
f  "A  ‘(f-v-y  w)  ,  1  “A  •  (.Z  -  y  v-y  w)  .  (2.11a,b) 

In  (2.9)  p  Is  the  reference  value  of  p  and  In  (2.11)  A^  are  the 

reciprocal  vectors  of  k.  Equation  (2.9a)  represents  conservation  of 

~1 

mass,  (2.9b,c)  represent  Che  balance  of  linear  momentum,  (2.9d,e) 
represent  Che  balance  of  director  momentum,  (2.9f,g)  represent  balances 
of  entropy,  and  (2.10a,b)  represent  balance  of  angular  momentum.  These 
equations  must  be  supplemented  by  an  energy  eo^'aclon  and  constitutive 
equations.  It  was  shown  in  [2]  that  the  energy  equation  for  a  linear 
chermoelastlc  shell  Is  satisfied  provided  that 


'o  ae 


9(1; 


ap 


a<P 

a(>  * 


M 


la 


Po  aic 


la 


(2.12a,b,c) 

(2.12d,e) 
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Additional  terms  appear  In  the  energy  equation,  but  they  are  of  higher 
order  and  thus  are  neglected  for  the  linear  theory. 


Now  we  confine  attention  to  a  thermoelastic  that  Is  Isotropic  In 
its  reference  configuration  and  specify  linear  constitutive  equations  of 
the  form 


2  a’®  f  «,(a“^  *  A-^  A^'')]e,p 

.  .  [,,  A»»  a''*  .  ..(A-"  A<^  *  a’^ 

*  "3  1'.  Ifp  *  “S  "Ba  ''3P  +  2  “9  T3 

-  2  ^  f3'«  -  «o>  -  2  h  -  2  h 

-  -  2  6  9^)  -  Sj  -  2  3,  9 

£  *  ■  s  •  El  “  ■  '*1  &1  ’ 

^  -  0  ,  Po  ?!  -  Po  ♦  , 


'^y6 


<l> 


(2.13a) 


(2.13b,c) 


(2.13d,e) 


where  the  coefficients  a^-a^,  Og,  a^,  ®q»  b^^,  b^  In  (2.13)  are 

constants  and  where  the  temperature  gradients  g  and  g  are  defined  by 

g  -  9  A®  ,  g,  -  (t.  A®  .  (2.14a,b) 

»  ,a  ~  *1  ,a£  ~ 

The  constitutive  assumption  In  (2.13a)  Is  slightly  different  from  that 
used  in  [2],  with  the  main  difference  being  that  here  tp  is  a  function  of 
instead  of  the  kinematic  variable  (see  additional  comments  in 
Section  5).  It  follows  from  (2.12)  and  (2.13)  that 


.®P  ^y6 


A®P  aT^°  +  a2(A®^  A^®  +  A®®  A^Me 


y6 


A«P/ 


+  atg  A-"  Y,  -  Pi  A'‘^(0  -  9  ) 


(2.15a) 


k“  -  «3  Yj  .  I"’  •  «4  Y3  “9  -  V  • 


-  [.,  a''^  *  aAa’’  A^  *  a“«  A^l')]. 


.«e 


«  ,Jic^g  -  ?2  ^  <*  ' 


u^“  a“P 

M  -  ap  A  <3p  . 


Po  “  pQ  Y3  +  Pi  a“^  e^p  +  P3(a  -  e^)  +  P5  . 

ap 


Po  ni  -  p2  A  ic^p  +  p^  <t> 


(2.15b,c) 

(2.15d) 

(2.15e) 

(2.15f) 

(2.15g) 


In  general,  these  constitutive  equations  oust  be  further  restricted  by 
statements  of  the  second  law  of  thermodynamics  [2].  For  the  thermo- 
elastic  shell  considered  here,  these  restrictions  reduce  to 


a>0  ,  b,  >0  ,  b,  >0  ,  p,  >0  .  (2 .16a,b,c,d) 

Before  closing  this  section,  we  recall  [1,2]  that  this  theory, 
which  Is  developed  by  direct  approach,  may  be  brought  Into  a  one-to-one 
correspondence  with  the  three-dimensional  theory  by  assuming  that  the 
position  vector  r*  of  a  point  In  the  shell  and  the  temperature  field  0* 
admit  the  representations 

r  -  r  (9®,9^,t)  -  r(9°',t)  +  9^  d(9®,t)  ,  (2.17a) 

9*  -  0*(9®,9^,t)  -  9(9®, t)  +  0^  0(0®, t)  ,  (2.17b) 

3 

where  9  Is  a  coordinate  through  the  thickness  of  the  shell.  For  a 
shell  of  constant  thickness  h,  we  may  choose  the  reference  surface  of 
the  shell  to  be  the  middle  surface  and  define  the  top  surface  SP"*"  of  the 
shell  by  9^  ■  h/2  and  the  bottom  surface  ap~  by  0^  -  -h/2.  With  reference 
to  the  three-dimensional  theory,  we  recall  [2]  the  definitions* 


*0ur  use  of  the  symbol  for  the  base  vector  should  not  be  confused 
with  the  use  of  the  same  symbol  for  the  temperature  gradient  In 
(2.Ub). 
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* 

11  1/2 

Si  - 1  >  Si  *  ’  ^1  *  «  *  [Si  S2  S3’  >  ^  '  (2.18a,b,c) 
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«lj  *  Si  *  &j  ’  ”  S^  *  S^  .  (2.18d,e 


where  and  are,  respectively,  the  base  vectors  and  their 
reciprocals  In  the  present  configuration  and  where  6^  Is  Che  Kronecker 
symbol.  We  also  recall  [2]  the  following  relationships  for  the  linear 


theory 


p*  0^/2  ^q3  ^ 

-h/2  ° 


Xy"  .  P*  de’ 


(N  -  1,2)  , 


-h/2 


\f  =-  /  p*  f*  d0^  +  b"*"  t'*’  +  b"  t"  , 

~  -h/2  ~  ~  • 


h/2 


U  -  /  p*  f*  0^  d0^  +  (|)  b"^  t"^  -  (y)  b‘  t"  , 


*1/2*  3  ++  -  - 

Xs  -  /  p^  s  dB-’  -  B  k  -  B  k  , 


-h/a 


Xs  -  /  ^  p*  G^^^  s*  9^  d9^  -  (y)  b"^  k"^  +  (y)  b"  k" 

^  -h/a  °  ^  ^ 


(2.19a) 


(2.19b) 


(2.19c) 


(2.19d) 


(2.9e) 


(2.19f) 


where  p  ,  f  ,  s  are,  respectively,  the  three-dimensional  mass  density 
(mass  per  unit  volume  In  the  reference  configuration),  specific  body 
force,  and  specific  rate  of  entropy  supply;  t^  and  t  are,  respectively, 
the  surface  tractions  applied  to  the  major  surface  5?"^  and  5P~;  k'*’  and 

k~  are,  respectively,  the  entropy  fluxes  applied  to  the  major  surfaces 

+  —  1 /7  1/2 

5P  and  5P  ;  and  G  '  Is  the  reference  value  of  g  .  It  follows  from 

(2.4),  (2.17),  and  (2.18)  that 
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(2.20) 


1/2  l/2r,  a  ^  ,-3.2  ,1  2  2  1., 

G  "A  [1  ~  0  +  (0  )  (Bj^  Bj  “  Bj^  ®2^J 


and  from  [2]  that  the  functions  B*^  and  B~  take  the  values 


.2  „1, 

>  ■  2  “o  ^  “2  ■  ^ 

Bl/2r,  .  h  .  IL  .2  «2 


B--A^'^[1+|b^+^(b{b^-B^  B^)]  . 


(2.21a) 

(2.21b) 


Also,  we  note  that  for  the  linear  theory  the  quantities  0^k  and  0  k 
represent,  respectively,  the  heat  fluxes  measured  positive  for  heat 
flowing  out  of  the  surfaces  and  5P~. 


3.  Determination  of  Constitutive  Coefficients 

In  this  section,  we  mainly  recall  results  that  were  obtained  In  [1] 
for  most  of  the  mechanical  constitutive  coefficients.  Confining 

* 

attention  to  a  shell  of  uniform  thickness  h  and  uniform  density  In 
Its  reference  configuration,  we  may  substitute  (2.20)  Into  equations 
(2.19a,b)  to  obtain  the  expressions 


^  *1/2  ,  *  .  .1/2. r-  ^  h^  ,_1 

X  -  A  -  (p^  h  A  )[!  +  12  ^®1 

A  1  /  *  1.  *l/2*/h^' 

Xy  -  -  (p^  h  A  )(^2  ®c,)  » 

^  2  ,  *  .  ,l/2.,h^.r,  ^  3  h^  ,_1  _ 
Xy  -  (p^  h  A  )('i2)[l  +  *^o“  ^®1  ® 

which  determine  the  reference  density  p^ 


2 

2  “ 
and 


-  bJ  B^)]  ,  (3.1a) 

(3.1b) 

bJ  B^)]  ,  (3.1c) 

the  Inertia  coefficients 


Values  for  most  of  the  constitutive  coefficients  were  obtained  In 
[1]  by  comparing  results  of  the  Cosserat  theory  with  exact  three- 
dimensional  results.  Even  though  the  constitutive  equations  (2.13)  and 
(2.15)  are  postulated  for  shells.  It  suffices  to  evaluate  most  of  the 
constitutive  coefficients  by  considering  solutions  of  plate  problems. 

By  solving  the  Isothermal  problem  of  a  plate  subjected  to  uniaxial 


C-12 


stress  (or  resultant  force)  and  that  for  a  plate  subjected  to  hydro¬ 
static  pressure  it  may  be  shown  [1]  chat 


a 


1 


v(l  -  V)  p 
(1  -  2v)  ^ 


“4 


(1  -  v)^ 

(1  -  2v) 


C 


ph 


(3.2a,b) 


(3.2c,d) 


where  E  Is  Young's  modulus,  v  is  Poisson's  ratio,  and  p  Is  one  of 
Lame's  constants,  all  associated  with  Che  three-dimensional  material. 
Similarly,  the  Isothermal  problem  of  pure  bending  of  a  plate  may  be 
solved  to  obtain  Che  results 


a,  -  vB  ,  a  -  B  ,  B- - — - z-  .  (3.3a,b,c) 

^  ^  2  12(1  -  v^) 

The  thermal  coefficients  and  the  remaining  mechanical  coefficients  will 
be  determined  In  the  next  two  sections. 


4.  Determination  of  the  Thermal  Coefficients 


To  determine  the  chermoelasClc  coefficients  ^21  ^  follow 
[1]  and  consider  the  problems  of  free  thermal  expansion  of  a  plate  and 
free  thermal  bending  to  obtain 


P 


o 


Pi 


* 

Eh  g 
1  -  2v 


3  * 

R  .  g 

^2  12(1  -  v) 


(4.1a,b) 


where  a  la  the  coefficient  of  linear  thermal  expansion  associated  with 
the  three-dimensional  material. 


The  remaining  thermal  coefficients  were  determined  in  [2]  by  direct 
Integration  of  the  constitutive  equations.  Here,  Instead,  we  determine 
these  coefficients  by  comparison  with  exact  solutions  of  a  rigid  heat 
conducting  plate.  Specifically,  we  consider  the  problem  of  a  plate  that 
Is  Initially  at  uniform  temperature  9^  and  that  Is  subjected  to  a 
uniform  heat  flux  on  Its  top  surface,  zero  heat  flux  on  Its  bottom 
surface,  and  no  external  supply  of  heat  (s*  ■  0).  The  quantity  q"*"  is 
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■  •  ..  •*  ^  *  .*  ■  .  ■ 


taken  to  be  positive  when  heat  flows  out  of  the  plate.  Using  k'*’  > 
and  k~  ■  0,  the  solution  of  equations  (2.9f,g)  may  be  written  in  the 
form 


9-0. 


This  Cosserat  solution  may  be  compared  with  the  exact  solution  recorded 

•k 

in  [4,  p.  112]  by  defining  the  average  temperature  6  and  the  average 

*  avg 

temperature  gradient  It*  direction  normal  to  the  plate's 

surface.  It  follows  that 


*  1^/2* 

0  =9+r/  (0-9  )d9-^  -  0 

avg  o  h  o  o 


-l^-V  . 

\P^  Ch/ 


(4.3a) 


exp  ( 

CD 

1  _  96  _ 

(2n-l)\n^t 

12  *  3  3  i 

5  -4-  /  (9  -  0^)0^  d0^  -  -  I 

h^  -h/2 

*  ^2 

Po  P  ^ 

i21l/ 

4 

L  X  n-l 

(2n  -  1)^ 

(4.3b) 


where  K  is  the  thermal  conductivity  and  c  is  the  specific  heat  at  zero 
strain.  For  most  practical  purposes,  we  may  retain  only  the  first  term 
in  the  series  in  (4.3b)  so  that  the  Cosserat  solution  and  the  exact 
solution  have  nearly  the  same  form  (note  that  96/ii^  -  0.9855).  Equating 
these  solutions,  we  obtain 


^3 


* 

p  c  h 

10_ 


P4 


* 

p  c  h 

9. 


(4.4a,b,c) 


Next,  we  compare  with  exact  steady  state  solutions  for  a  rigid  heat 
conductor  of  the  forms  9*  -  a9^  and  9*  -  a0^9^,  where  a  is  a  constant. 
This  comparison  yields 


h  - 


Kh' 


12  9 


(4.5a,b) 
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The  last  thermal  coefficient  corresponds  to  the  arbitrary  constant 
reference  value  of  the  entropy  and  therefore  cannot  be  specified. 


Apart  from  minor  differences  In  sign,  which  are  caused  by  our 
constitutive  assumption  (2.11a),  all  of  the  values  for  the  thermal 
coefficients  determined  here  are  the  same  as  those  proposed  In  [2], 
except  that  for  0^.  The  value  for  0^  proposed  In  [2], 


P4 


* 

P  h 

12  e 


(4.6) 


was  obtained  by  the  method  of  direct  Integration  of  the  constitutive 
equations  Instead  of  by  the  method  of  comparison  with  exact  solutions. 
Because  we  ultimately  require  the  Cosserat  theory  to  reproduce  exact 
results  with  as  little  error  as  possible,  we  adopt  the  latter  approach 
and  specify  0^  by  (4.4c)  Instead  of  (4.6). 

5 .  Determination  of  the  Mechanical  Coefficients  and  gR 

The  discussion  of  constitutive  equations  In  [1]  emphasized  that  In 
general,  the  constitutive  coefficients  are  not  constants.  In  other 
words,  values  that  are  obtained  by  comparing  Cosserat  solutions  with 
exact  solutions  of  one  class  of  problems  may  not  be  the  same  as  values 
obtained  by  considering  another  class  of  problems.  For  example,  values 
that  predict  accurate  results  for  certain  quantities  In  static  problems 
may  be  different  from  those  that  predict  accurate  results  In  dynamic 
problems.  Nevertheless,  from  a  practical  point  of  view  we  need  to 
specify  values  for  the  constitutive  coefficients.  By  way  of  background, 
we  recall  [1]  that  the  coefficients  and  Cg  associated  with  tangential 
shear  deformation*  have  been  specified  by 


*For  a  shell  (B^g  ♦  0)  the  coefficient  ag  relates  in  (2.15e)  to 
tangential  shear  (see  2.61). 


5  . 

03  -  g  Hh 


7  v3 

*8  “  l20 


(5.1a,b,c) 


2 

or  ®3  “  72 


> 


The  values  (5.1a,c)  can  be  motivated  by  assuming  a  form  for  the  stress 
distribution  through  the  thickness  of  a  plate  and  developing  an  approx¬ 
imate  expression  for  the  strain  energy  function  by  Integration  [1,5]. 
However,  the  value  (5.1b)  may  be  obtained  by  comparing  with  a  solution 
for  a  vibrating  plate  [8].  Thus,  we  would  expect  that  either  of  the 
values  (5.1a,b)  for  03  (those  values  are  very  close  to  each  other)  would 
be  appropriate  for  dynamical  problems.  However,  as  was  pointed  out  In 
[1],  the  value  (5.1c)  for  Og  has  not  been  validated  by  comparison  with 
any  exact  solution.  Here,  we  determine  different  values  for  and 
by  comparing  with  exact  solutions  of  two  static  problems. 

First,  we  determine  by  considering  the  static  Isothermal  problem 
of  simple  shear  In  the  6  -  9  plane  of  a  plate  In  the  absence  of  body 

force  (Jt*  ■  0).  For  this  problem,  we  specify  t^  -  r  ,  £  ■  -  t  , 
where  t  Is  the  shear  stress  applied  to  the  major  surfaces  of  the  plate. 
It  follows  that  the  Cosserat  solution  Is  an  exact  solution  If  we  specify 

O3  -  ph  .  (5.2) 

A  discussion  of  this  result  will  be  given  at  the  end  of  this  section. 

Next,  the  coefficient  Ug  will  be  determined  by  considering  the 
static.  Isothermal  problem  of  a  circular  cylindrical  shell  of  radius  R 
and  thickness  h  with  Its  Inner  surface  held  fixed  and  Its  outer  surface 
rotated  by  an  amount  xh  (Figure  C.l).  Let  be  an  orthonormal 
coordinate  system  with  parallel  to  the  circumferential  tangent  to  the 
shell  and  (gg  parallel  to  the  outward  normal  to  the  shell  (Figure  C.l). 
Further,  let  r*  be  the  distance  of  a  point  from  the  symmetry  axes  of  the 
shell,  be  the  displacement  In  the  jg^  direction,  and  t^g  be  the  (1,3) 
component  of  the  Cauchy  stress  tj^j.  Then  the  exact  three-dimensional 
solution  yields 

r*  -  R  +  9^  ,  (5.3a) 
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FIGURE  C.1  CROSS  SECTION  OF  A  CIRCULAR  CYLINDRICAL  SHELL 
WITH  THE  INNER  SURFACE  9P'  HELD  FIXED  AND  THE 
OUTER  SURFACE  dP*  ROTATED  BY  AN  AMOUNT  Kh 

For  this  solution  R  =  h. 


(5.3b) 


l/i  .  h  1  r  *  h \2  R  i 

+  2R^  ^  "  2R^ 


tl3(t‘)  -  t„(r*)  -  1.(1  .^)(1  -fe) 


"31' 


2R^ 


2«'  (^.,2-  • 


(5.1c) 


where  all  other  components  of  the  displacement  and  stress  *^13  vanish. 
Thus,  the  displacements  u,^  and  ju  ,  and  traction  vectors  and  £*  on  the 
top  surface  and  bottom  surface  bP”,  respectively,  are  given  by 


+ 

u 


t 


•=13^®  2^£i  “ 

-  -  j)®!  - 


£i  »  H  ■  T^~i  "  ^  ’ 

(5.4a,b) 

fi  - 

^  2R^  .+ 

.  ic  e,  «  t  e,  , 

(1  -  If)  "  " 

(5. Ac) 

+  -  -  t‘ .1  , 

(5. Ad) 

where  t"*"  and  t”  are  the  values  of  the  stress  ti3  on  the  top  and  bottom 
surfaces,  respectively.  The  solution  u*(r*)  Is  plotted  In  Figure  C.l 
for  the  thick-shell  case  where  R  •  h.  Notice  that  since  u*  Is  nearly 
linear  In  r  ,  we  should  expect  the  Cosserat  solution  to  predict  accurate 
results  even  In  the  limit  of  a  thick  shell. 

Now,  the  geometry  of  a  cylindrical  shell  of  radius  R  is 
characterized  by 


A  „ 

~a,3 


»  0  ,(5.5a,b,c,d) 


all  other  B  „  ■  0  , 


all  other  B^  ■  0  , 


(5.5e,f) 

(5.5g,h) 


where  r^o  Is  the  Chrlstoffel  symbol.  Using  (2.5),  (2.17a)  and  (5.5),  we 

®p 

realize  that  the  conditions  (5.Aa,b)  yield 
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.h,.  ~  h,_ 

u  -  u+  -:r6“iche,  ,  u  -  u-  ttO-O 


Consistent  with  (5.6)  we  may  take 


(5.6a,b) 


1  hic  2  3  .  .1  ic  ,2  ,3  -  -v 

u--^,  u-u-0  ,  6  6  “6  -  0  ,  (5.7a,b,c,d) 


so  that 


Yl  -  Rk(1  -  . 


(5.8a,b) 


'.S  ■  “  •  ■  ''3  ■  "  •  'as  ■  “  •  '32  ■  "  •  <5-8<=.d...() 

It  follows  that  the  only  nontrivial  equations  of  equilibrium  become 


P  (f^  +  7  -  0  ,  0  -  p  ^ 


(5.9a,b) 


where 


(5.10a) 


(5.10b) 


k‘  .  -  5^)  ,  „31  .  .  .  1-)  .  (s.lOc.d) 

Substituting  (5.10)  Into  (5.9)  and  comparing  the  result  with  (5.4),  we 
realize  that  the  Cosserat  equations  will  predict  exact  values  for  t"^  and 
t“  If  and  Og  satisfy  the  equation 


a,  +  -r  ■  ph  . 

Now  If  we  adopt  the  specification  (5.2),  we  deduce  that 


(5.11) 


-  0  . 


(5.12) 


Ue  are  now  In  a  position  to  discuss  the  significance  of  the  results 

(5.2)  and  (5.12).  Since  the  value  of  given  by  (5.2)  was  obtained  by 

comparing  with  the  exact  solution  of  the  static  simple  shear  problem.  It 

follows  that  If  were  specified  by  (5.1a,b)  Instead  of  (5.2),  the 

Cosserat  equations  would  necessarily  predict  Incorrect  results  for 

simple  shear  which  would  be  undesirable.  However,  one  could  question 

whether  It  Is  of  prime  Importance  for  a  shell  theory  to  predict  results 

of  simple  shear — which  requires  surface  tractions  to  be  specified  on  the 

major  surfaces — when  many  applications  of  shell  theory  consider  shells 

with  free  major  surfaces.  In  response  to  this  question,  we  note  that 

two  Important  problems  of  a  plate  iflth  free  major  surfaces  have  been 

solved  for  arbitrary  values  of  a^.  Problem  A  [1,  Sec.  24]  considers 

pure  twist  of  a  plate  and  problem  B  [7]  considers  pure  bending  of  a 

plate  with  a  circular  hole.  If  the  value  [5.2]  for  a-^  Is  adopted 

Instead  of  (5.1a),  then  It  can  be  shown  that  the  solution  of  problem  A 

Is  slightly  Improved  and  the  solution  of  problem  B  Is  only  slightly 

1  /2 

modified  with  q  being  replaced  by  (6/5)  '  q  ■  1.10  n  In  the  formulas  In 

[7]. 

It  Is  also  of  Interest  to  note  that  If  we  were  to  specify  alterna¬ 
tive  constitutive  equations  by  replacing  In  (2.13a)  with  then 

on 

for  the  example  associated  with  Figure  C.l,  the  quantity  M-"-  would  vanish 

and  equations  (5.9)  would  again  yield  the  result  (5.2).  With  reference 

to  this  same  example,  we  observe  that  there  Is  an  Inconsistency  between 

the  specifications  (5.1)  and  the  result  (5.11).  In  particular.  It  Is 

clear  that  If  either  (5.1a)  or  (5.1b)  are  substituted  Into  (5.11),  we 

2 

conclude  that  ag  Is  proportional  to  phR  ,  which  Is  significantly 
different  than  (5.1c). 

In  conclusion;  it  appears  that  for  static  problems  It  Is  better  to 
specify  ttj  by  (5.2)  Instead  of  (5.1a,b)  and  by  (5.12)  Instead  of 

(5.1c).  This  Is  because  with  this  specification,  the  Cosserat  theory 
predicts  accurate  results  for  all  four  static  problems  considered  In 
this  section  without  the  Inconsistency  described  above.  However,  If 
comparison  with  the  dynamic  problem  considered  In  [6]  Is  of  prime 


Importance,  then  It  is  best  to  specify  by  (5 •lb)  Instead  of  (5.2). 
In  conclusion,  we  recall  the  observation  In  [1]  that  the  constitutive 
coefficients  for  shell  theory  are  not  constants  In  the  sense  that  the 
best  values  for  these  coefficients  may  be  problem  dependent. 
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Appendix  D 

HEAT  CONDUCTION  IN  PLATES  AND  SHELLS  WITH 
EMPHASIS  ON  A  CONICAL  SHELL 


(Submitted  for  publication  In  the  International 
Journal  of  Solids  and  Structures) 


ABSTRACT 


This  paper  Is  concerned  with  analyzing  heat  conduction  In  rigid 
shell-llhe  bodies.  The  thermal  equations  of  the  theory  of  a  Cosserat 
surface  are  used  to  calculate  the  average  (through-the-thlckness) 
temperature  and  temperature  gradient  directly.  Specific  attention  Is 
focused  on  a  conical  shell.  The  conical  shell  Is  particularly  Interest 
Ing  because  It  has  a  converging  geometry  such  that  the  shell  near  Its 
tip  Is  "thick”  even  though  the  shell  near  Its  base  may  be  "thin.” 
Generalized  constitutive  equations,  which  Include  certain  geometrical 
features  of  shells,  are  developed  here  In  a  consistent  manner.  These 
equations  are  tested  by  considering  a  number  of  problems  of  plates, 
circular  cylindrical  shells,  and  spherical  shells,  and  comparing  the 
results  with  exact  solutions.  In  all  cases,  satisfactory  results  are 
predicted  even  In  the  thlck~shell  limit.  Finally,  a  problem  of 
transient  heat  conduction  In  a  conical  shell  Is  solved.  It  Is  shown 
that  the  thermal  bending  moment  produced  by  the  average  temperature 
gradient  Is  quite  severe  near  the  tip  and  It  attains  Its  maximum  value 
In  a  relatively  short  time. 
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Introduction 


Most  aerospace  structures  are  compositions  of  structural  components 
that  can  be  modeled  as  shell-llke  bodies.  For  various  reasons.  It  Is 
desirable  to  determine  the  thermomechanical  response  of  these  shell-llke 
bodies  to  thermal  and  mechanical  loads.  Within  the  context  of  classical 
linear  shell  theory,  the  temperature  distribution  Influences  the 
mechanical  response  of  the  shell  through  the  resultant  thermal  force  and 
resultant  thermal  moment.  For  an  elastic  shell,  the  thermal  force  Is 
related  to  the  average  ( through- the- thickness)  temperature  and  the 
thermal  moment  Is  related  to  the  average  temperature  gradient  by 
constitutive  equations. 

If  the  strain  rate  In  a  given  problem  Is  sufficiently  small,  then 
the  thermal  and  mechanical  problems  uncouple  In  the  sense  that  the 
temperature  field  may  be  determined  by  solving  equations  for  a  rigid 
heat  conductor.  Then  the  resulting  temperature  field  may  be  used  to 
calculate  the  thermal  force  and  thermal  moment,  which  provide  thermal 
loading  for  determining  shell  deformation. 

In  this  paper,  we  confine  attention  to  determining  the  temperature 
distribution  In  a  shell-llke  body  that  Is  treated  as  a  rigid  heat 
conductor.  Although  the  temperature  distribution  can  be  determined  by 
attempting  to  solve  the  three-dimensional  heat  conduction  equation,  this 
approach  has  two  major  disadvantages.  First,  because  the  thermal  loads 
for  shell  theory  depend  only  on  the  average  temperature  and  temperature 
gradient,  much  of  the  details  calculated  by  this  approach  are  not  of 
prime  Importance.  Second,  because  the  heat  conduction  equation  admits 
separable  solutions  for  only  a  limited  number  of  geometries.  It  Is 
exceedingly  difficult  to  obtain  analytical  solutions  for  many  typical 
shell  geometries.  This  latter  problem  has  been  addressed  In  [1],  where 
equations  are  developed  to  calculate  an  approximate  temperature 
distribution  In  shells  of  revolution. 

We  take  a  different  approach  and  use  thermal  equations  for  shells 
that  have  recently  been  developed  [2]  to  predict  the  average  temperature 
and  temperature  gradient  directly.  These  equations  are  based  on 


modeling  the  shell-llke  body  as  a  Cosserat  surface.  Details  of  this 
theory  may  be  found  In  [2,3].  Specifically,  the  objective  of  this  paper 
Is  to  determine  the  average  temperature  and  temperature  gradient  In  a 
conical  shell  (Figure  D.l),  which  Is  a  basic  aerospace  structure.  The 
conical  shell  Is  particularly  Interesting  because  It  has  a  converging 
geometry  so  tha  the  shell  near  Its  tip  Is  necessarily  “thick"  even 
though  the  shell  near  Its  base  may  be  "thin."  For  this  reason.  It  Is 
questionable  whether  any  shell  theory  can  accurately  predict  results  for 
the  critical  tip  region.  Here  It  Is  shown  that  with  appropriate 
constitutive  equations,  the  Cosserat  theory  Includes  enough  of  the 
geometry  of  the  shell  to  predict  relatively  accurate  results  for  the 
conical  shell. 

It  Is  not  a  trivial  matter  to  develop  equations  for  shells  that 
produce  reasonable  results  In  the  thick-shell  limit.  For  example,  we 
recall  that  the  equations  In  [1]  were  developed  by  writing  the  heat 
conduction  equation  In  a  form  appropriate  for  shells  and  then  neglecting 
quantities  multiplied  by  higher  powers  of  the  ratio  of  the  thickness  to 
radius  of  curvature.  Even  though  these  equations  are  more  complicated 
than  the  Cosserat  equations  In  that  details  of  the  through-the-thlckness 
temperature  distribution  are  calculated,  too  much  of  the  shell  geometry 
has  been  neglected,  and  hence  they  predict  Inaccurate  results  In  the 
thick-shell  limit.  The  predictions  of  the  equations  In  [1]  are  compared 
with  the  more  accurate  predictions  of  the  Cosserat  theory  for  the  thick- 
shell  problems  considered  In  Sections  4  and  5. 

In  the  following  sections,  we  discuss  the  basic  equations  of  the 
Cosserat  theory  and  then  solve  a  number  of  problems.  To  develop 
confidence  In  the  predictions  of  the  Cosserat  theory  In  the  base  region 
of  the  conical  shell,  we  solve  various  problems  for  a  plate  and  compare 
with  exact  solutions  In  [A].  These  problems  examine  the  effect  of  the 
three  types  of  boundary  conditions  on  the  major  surfaces  of  the  plate: 
specified  heat  flux,  specified  temperature,  and  radiation.  Next,  to 
develop  confidence  In  the  predictions  of  the  theory  In  the  tip  region  of 
the  conical  shell,  we  use  the  same  equations  to  solve  specific  problems 
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for  a  solid  circular  cylinder  and  a  solid  sphere,  and  compare  the 
results  with  exact  solutions.  Finally,  after  having  developed 
confidence  In  the  predictions  of  the  theory  In  both  the  tip  and  base 
regions  of  the  conical  shell,  we  solve  a  specific  heat  conduction 
problem  for  a  conical  shell. 


2 .  Basic  Equations 

Let  the  material  points  of  the  Cosserat  surface  C  be  Identified  by 
means  of  a  system  of  convected  coordinates  9‘’‘(a  >  1,2}  and  let  the  two- 
dimensional  region  of  space  occupied  by  the  material  surface  In  the 
present  configuration  at  time  t  be  denoted  by  c.  Further,  let  the 
vector  valued  function  x  define  the  position  of  a  material  point  of  the 
surface  C  and  at  each  such  point  define  the  vector  valued  function  i, 
called  the  director,  and  the  two  temperature  fields  3  and  each 
referred  to  the  present  configuration.  Then  a  therraomechanlcal  process 
of  the  Cosserat  surface  Is  defined  by 

r  -  r(9®,t)  ,  d  -  d(0“,t)  ,  [a^  *2  °  »  (2.1a,b,c) 

9  -  9(9®, t)  ,  (9  >  0)  ,  4.  -  ♦(9®,t)  ,  (2.1d,e,f) 


where  the  tangent  vectors  ^  and  the  unit  normal  vector  ^3  are  defined 
by 


a 

~a 


ar 


99 


~a 


•  a 


"3 


,  a,  •  a-  -  1 
’  '^3  ~3 
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(2.2a,b,c,d) 


and  the  condition  (2.1c)  ensures  that  the  director  Is  nowhere  tangent  to 
c.  Also,  In  the  above,  9  represents  the  average  (through-the-thlckness) 
temperature  In  the  shell  and  ^  represents  the  average  temperature 
gradient. 

In  the  reference  configuration,  we  assume  that  the  shell  has 
uniform  thickness  h  and  Is  at  uniform  temperature  9^.  Then,  the 
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.  '•j. 

■>  'v- '  ■ '  -  V  ■  -  ■  -  '  .'V  w'-  ."'V" 

r.'"’-'*'--'' v.-'v'v'  v-'-  ' 


reference  values  of  the  various  kinematic  quantities  may  be  denoted  by 


r»R  ,  d^D^Aj  t  a^^Aj  ,  ,  (2.3a,b,c,d) 

9-0  ,  4)  -  0  ,  (2.3e,f) 

o 

where  J,  and  A^^^  depend  on  the  coordinate  9®  only.  For  the  rigid 
heat  conductor  considered  here,  there  is  no  distinction  between  the 
reference  and  the  present  configurations;  hence  equations  (2.3a-d)  hold 
for  all  time.  Further,  all  tensor  quantities  will  be  referred  to  the 
base  vectors  their  reciprocals  defined  by 

A,  •  A^  -  ,  (2.4) 

where  5^^  Is  the  Kronecker  symbol. 

Let  P,  bounded  by  the  closed  curve  3P,  denote  the  region  occupied 
by  an  arbitrary  material  portion  of  the  surface  c  and  let  be  the  unit 
outward  normal  to  dP.  Using  the  notation  of  [2],  we  define  the  follow¬ 
ing  quantities:  the  positive  mass  density  (mass  per  unit  area  of  P)  In 
the  reference  configuration  -  Pq(9‘’^);  the  specific  (per  unit  mass  of 
P)  entropies  q  -  q(9®,t)  and  -  qi(9®,t);  the  specific  Internal  rates 
of  production  of  entropy  ^  -  5(9®, t),  **  Ci(9®»t),  and  -  5j^(9*,t); 

the  entropy  fluxes  k  -  k(9®,t;  v)  and  kj^  -  k^(0“,t;  v),  each  per  unit 
length  of  the  curve  9P;  the  specific  external  rates  of  supply  of  entropy 
8  -  8(9®, t)  and  s^  -  8^(0®, t);  the  specific  Internal  energy  e  -  e(0®,t); 
and  the  specific  Helmholtz  free  energy  (|;  -  (|;(9®,t)  = 
e  -  9q  -  ^q|^.  With  suitable  continuity  assumptions.  It  can  be  shown 
that  [2,3] 

k  -  £  •  V  -  p“  ^  ”  Pi  ''a  ’  (2.5a,b) 


^Throughout  the  text  Greek  Indices  have  a  range  (1,2)  and  Latin  Indices 
have  a  range  (1,2,3). 


o 


where  v  ■  A  •  v  are  the  components  of  the  normal  vector  v  and  where  we 
a  ~o  ~ 


use  the  usual  summation  convention  over  repeated  Indices*  Further,  with 
reference  to  the  energy  equation,  the  specific  external  rates  of  heat 


supply  r  and  r^;  and  the  heat  flux  vectors  £  and  are  defined  by 


r-0a  .  r^-((,Sj  ,  a-9£  . 


(2.6a,b,c,d) 


Now,  the  local  forms  of  the  balances  of  entropy  may  be  recorded  as 


pQ  ^  ”  Po^®  +  “  p“la  •  Po  “  Po^®l  ^1^  ■  Pl“l(x  *  (2.7a,b) 


where  a  dot  denotes  material  time  differentiation  and  where  a  bar 


denotes  covariant  differentiation  with  respect  to  the  metric  of  the 
shell  surface.  For  later  convenience,  we  recall  [2,3]  definitions  for 
the  metric  tensor  A-g,  and  Its  reciprocal  A'*^,  the  curvature  tensor  B_g, 


the  Christoff el  S3rmbol  r°p,  and  covariant  differentiation  In  the  forms 


*afl  ■  4.  •  ip  •  ■  4*  •  4^  .  ■  *..p  •  43  .  (2.8..b,  =  ) 


''ap  ■  i..P  •  i“  '  •  "“ip  ■  '>°.p  -  '■^p  >  (Z-S'i.'.f) 


where  a  comma  denotes  partial  differentiation  with  respect  to  6°^. 


Equations  (2.7)  must  be  supplemented  by  an  energy  equation  and 
constitutive  equations.  It  follows  from  [2]  that  the  energy  equation 
for  a  rigid  thermoelastic  shell  Is  satisfied  provided  that 


9(ti 

■  d(ti  *■ 


(2.9a,b) 


Pq  0  5  +  Po  0  ^1  +  P  •  &  +  £l  •  ll  -  0 


(2.10) 


where  the  temperature  gradients  g  and  are  defined  by 


&  -  ®.a  *  Si  -  •  (2.11a.b) 

Confining  attention  to  a  rigid  shell  that  Is  thermally  Isotropic,  we 
specify  constitutive  equations  In  the  form 

2  (1,  -  -  ^3(9^  -  2  0  6^)  -  9^  0^  -  2  P3  0  ,  (2.12a) 

P  "  -  a  g  ,  p  -  -  b  g  ,  (2.12b,c) 

~  o  ~  ~l  1  ~l 

pQ  ®  ^  “  *0  S  *  &  ’’l  Si  *  Si  +  ^2  '  (2.12d) 

Pq  ”  Pq  ^1  ■  ■  **2  (2.12e) 

where  a^,  b^,  b2,  ^3  -  Pj  are  constants.  Substituting  (2.12a)  Into 
(2.9),  we  have 

pQ  T)  ■  PjCQ  -  +  P5  .  Pp  1^1  "  P4  4'  •  (2.13a,b) 

The  form  of  the  constitutive  equations  (2.12)  represents  a  slight 
generalization  of  those  Introduced  In  [2]  for  the  linear  theory.*  These 
equations  are  chosen  to  automatically  satisfy  the  reduced  energy 
equation  (2.10)  without  approximation. 

In  postulating  the  form  of  the  constitutive  equations  (2.12),  we 
tacitly  assume  that  constitutive  equations  that  are  valid  for  a  plate 
are  also  valid  for  a  shell.  In  the  discussion  In  Section  5,  we  observe 
that  certain  geometrical  features  of  the  shell  must  be  Included  In  the 
constitutive  equations  to  predict  relatively  accurate  results  for  a 
solid  sphere.  These  geometrical  features  of  the  shell  may  be  Introduced 


*The  sign  of  the  constants  P3,  pA  used  here  Is  opposite  that  used  In 

[2]. 


It 


O 
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by  appropriately  modifying  the  constitutive  equations  to  take  the  forms 
2  p*  h  ♦  -  -  -  2  e  9^)  -  pj  -  2  6  ,  (2.14i.) 


£--•„£  •  El  -  -  »! 


Po  ®  5  ^  •  ft  +  ^  &1  *  &1  +  (-^) 


Po  *'  ^1  “  Po  ^  ^1  “  ■  ^2  ♦  * 


b2  ♦  . 


(2.14b,c) 


(2.14d) 


(2.14e) 


P  P 

Pq  h  -  [PjO  -  Qq)  +  Pjl  ,  pQ  ’ll  “  P4  *  »  (2.14f,g) 

Po  ”  Po 

* 

where  Is  the  constant  three-dimensional  mass  density  (mass  per  unit 
volume)  of  the  material  and  h  Is  the  constant  thickness  of  the  shell. 

The  constitutive  equations  (2.14f,g)  depend  on  the  geometry  of  the  shell 

* 

through  the  ratio  p^/p^  h  [see  equation  (2.16a)]. 

Within  the  context  of  the  general  theory,  the  constitutive 
equations  must  be  further  restricted  by  statements  of  the  second  law  of 
thermodynamics  [2].  For  either  of  the  sets  of  constitutive  equations 
(2.12)  and  (2.13)  or  (2.14),  these  restrictions  reduce  to 


«o  -  ”  ’  •  ^2  -  °  ’  P3  ^  ^  • 


(2 .15a,b,c,d) 


To  linearize  the  equations  presented  above,  we  assume  that  the  tempera¬ 
tures  (9  -  6g)  and  and  their  space  and  time  derivatives  are  small 
enough  that  quadratic  expressions  In  these  quantities  may  be  neglected 
relative  to  linear  expressions.  It  follows  from  (2.12d)  or  (2.14d)  that 
^  Is  of  higher  order  so  that  ^  may  be  set  equal  to  zero  In  (2.7a). 

Now,  we  recall  [2,3]  that  the  Cosserat  theory  developed  by  direct 

approach  may  be  brought  Into  a  one-to-one  correspondence  with  the  three- 

* 

dimensional  theory  by  assuming  that  the  position  vector  r  of  a  point  In 


/  A 


the  shell  and  the  temperature  field  9  admit  the  representations 

r*  -  r*(9“,0^,t)  -  £(0“,t)  +  0^  d(0“,t)  ,  (2.15a) 

9*  -  9*(0“,0^,t)  -  9(0“,t)  +  9^  ^(e“,t)  ,  (2.15b) 


where  0^  Is  a  coordinate  through  the  thickness  of  the  shell.  For  a 
shell  of  constant  thickness  h,  we  may  choose  the  reference  surface  of 
the  shell  to  be  the  middle  surface  and  define  the  top  surface  dP^  of  the 

shell  by  0^  ■  h/2  and  the  bottom  surface  &P~  by  0^  ■  -  h/2.  If  the 

* 

three-dimensional  mass  density  of  the  shell  Is  constant,  then  It  may 
be  shorn  that  [2,3,5]: 


-  Po  p*  q1/2  ^  (bJ  -  bJ  B^)]  , 

(2.16a) 

\s  -  \8  -  b"*"  k'*’  -  b"  k"  ,  Xs  -  /  p*  s*  d0^  ,  (2.16b,c) 

-h/2  ° 

Xs.  -  Xs,  -  (^]  B^  k"^  +  (^)  b‘  k"  ,  Xs  -  /  p*  s*  9^  d0^  , 

112  2  1  ° 

(2.16d,e) 


^1/2  l/2r,  3  a  ^  ,„3.2  .1  2  2  1.-, 

G  "A  [1-0  B^  +  (9  )  (Bj^  B^  ~ 


B*  -  -  J  11°  +  J-  (»} 


9  «9+y(>  ,  0  ”9-^ 


bJ  B^)J  , 


bJ  bJ)]  , 


(2.16f) 


(2.16g) 


(2.16h) 


(2.161,j) 


»ss 


.v.v--. 


where  a*  Is  the  three-dimensional  rate  of  entropy  supply;  k"*"  and  k~  are, 
respectively,  the  entropy  fluxes  applied  to  the  major  surfaces  dP^  and 
dP~;  Bq  are  the  mixed  components  of  the  curvature  tensor;  and  9^  and  6~ 

P  .  _ 

are,  respectively,  the  temperatures  on  the  major  surfaces  and  dP  . 

Also,  we  note  that  for  the  linear  theory 


0k^»q^,0k--q  ,  (2.17a,b) 

o  o 

where  q'*'  Is  the  heat  flux  measured  positive  for  heat  flowing  out  of  the 
surface  dP^  and  q~  Is  Che  heat  flux  measured  positive  for  heat  flowing 
Into  Che  surface  dP~. 

Most  of  the  constitutive  coefficients  were  evaluated  In  [2]  by 
direct  Integration  of  Che  three-dimensional  constitutive  equations.  An 
alternative  approach  was  taken  In  [5],  where  Che  coefficients  were 
evaluated  by  comparing  Cosserat  solutions  with  exact  three-dimensional 
solutions.  Except  for  Che  value  of  ^3,  the  results  In  [2]  and  [5]  are 
Che  same.  Here,  we  adopt  Che  results  In  [5]  and  specify 


,  Kh^ 

-12-5; 


* 


.  P/. 


p  ch 
9 


(2.18a,b,c) 


(2.18d,e) 


where  K  Is  the  thermal  conductivity  and  c  Is  the  specific  heat  at 
constant  strain  of  the  material.  The  coefficient  85  corresponds  to  the 
arbitrary  constant  reference  value  of  the  entropy  and  therefore  cannot 
be  specified.  Since  the  material  constants  K  and  c  are  positive,  we 
realize  from  (2.18)  chat  Che  restrictions  (2.15)  are  satisfied. 

Finally,  we  use  (2.14),  (2.16)-(2.18)  to  write  the  linearized 
version  of  equations  (2.7)  In  Che  form 


•  -1/2  +  +  -  -  2 
pQ  c  0  “  Pq  00  s  -  A  [B  q  -  B  q  ]  +  Kh  V  9 


(2.19a) 


D-12 


o 


•  ^  .•1/2  fh\  -  +  +  —  — - 

2 —  *•>  *  pQ  ®o  ®i  "  ^  (yj  I®  q  +  ®  q  ] 

n 

-  Kh  (|»  +  ^  7^  d)  ,  (2.19b) 

Po 

2 

where  the  Laplaclan  operator  V  3  is  defined  by 

7^  9  -  A®^  0.  o  -  a“^(0  a  -  r®.  0  )  .  (2.20) 

lap  ,ap  ap  ,a' 

3.  Plates 

In  this  section,  we  examine  the  validity  of  the  Cosserat  theory  in 
the  thin-shell  limit  by  considering  three  problems  of  heat  conduction  in 
a  plate.  These  problems  are  chosen  to  examine  the  effects  of  specifying 
heat  flux,  temperature,  or  radiation-type  boundary  conditions  on  the 
major  surfaces  of  the  plate.  For  each  of  these  problems,  we  neglect 
external  entropy  supply  (or  external  heat  supply)  and  consider  tempera¬ 
ture  fields  that  are  functions  of  time  only  so  that 

A  A 

s-0  ,  s^-0  ,  0  -  9(t)  ,  -  (Kt)  .  (3.1a,b,c,d) 

Further,  the  curvature  tensor  for  a  plate  vanishes.  Hence,  from 
(2.16)  we  deduce  Chat 

*  +  1/2  -  1/2 
•  P«  "  P«  ^  »  ®  "  ^  •  (3.2a,b,c,d) 

ap  o  o 

and  chat  equations  (2.19)  reduce  to 

*  .  +- 
ch  0  -  -  q  +  q  , 


(3.3a) 


Problem  1;  For  this  problem,  the  heat  flux  q"*"  is  specified  to  be 
constant  on  the  top  surface,  the  bottom  surface  is  insulated,  and  the 
plate  is  initially  at  uniform  temperature  3^.  Mathematically,  these 
conditions  are  characterized  by 


+ 

q  *  constant 


q  ■  0 


(3.4a,b) 


6-9  ,  ^-0  at  t-0  .  (3.4c,d) 

o 

Since  the  solution  of  equations  (3.3)  %d.th  the  conditions  (3.4)  was 
developed  in  [5],  «ie  merely  record  the  solution  in  the  nondlmensional 
form 


-  K(0  -  0  ) 
o 


h  q 


-  y  [1  -  exp(-  t)1 

q 


(3.5a,b) 


where  x  is  the  nondlmensional  time  parameter  defined  by 


t 


Kt 

■*  2  * 

Po  ch 


(3.6) 


Recall  from  [5]  that  the  constitutive  coefficients  were  chosen  by 
requiring  the  Cosserat  solution  to  compare  very  well  with  the  exact 
solution  recorded  in  [4,  p.  112]. 

To  exhibit  this  comparison  graphically,  we  have  used  (2.15b)  to 
plot  in  Figure  D.2  the  Cosserat  solution  (3.5)  together  with  the  exact 
solution  for  various  values  of  the  time  parameter  t.  The  dashed  lines 
in  Figure  0.2  have  been  taken  directly  from  [4,  Figure  15,  p.  113]  and 

X  -  0^  +  j  (3.7) 

so  that  X  •  0  locates  the  bottom  surface  9P~  and  x  -  h  locates  the  top 
surface  dP"^. 
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FIGURE  D.2  NORMALIZED  TEMPERATURE  IN  A  PLATE  OF  THICKNESS  h, 
WITH  ZERO  HEAT  FLUX  AT  x  =  0,  CONSTANT  HEAT  FLUX  q* 
(OUT  OF  THE  PLATE)  AT  x  =  h,  AND  UNIFORM  INITIAL 
TEMPERATURES 

The  numbers  on  the  curves  are  values  of  t  “  Kt/p*  ch^.  The  dashed 
lines  are  the  exact  solution  and  the  solid  lines  are  the  Cosserat 
solution. 
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Problem  2;  For  this  problem,  the  temperature  O’*"  Is  specified  as  0^  on 
the  top  surface,  the  heat  flux  q~  is  specified  to  be  constant  on  the 
bottom  surface,  and  the  plate  Is  Initially  at  uniform  temperature  @g> 
Mathematically,  these  conditions  are  characterized  by 


constant 


0-0 


and  ^ 


at 


With  the  help  of  (2.161)  condition  (3.8a)  yields 


(3.8a,b) 


(3.8c,d) 


«  -  ®o  -  2 


(3.9) 


It  is  Important  to  observe  here  that  by  specifying  0'*',  we  tacitly 
specify  0  In  terms  of  ^  through  equation  (2.161).  It  follows  that  it  Is 
not  possible  to  specify  Independent  Initial  values  for  0  and  such  as 
(3.8c,d).  In  other  words,  when  temperature  Is  specified  on  one  or  both 
of  the  major  surfaces,  we  must.  In  general,  modify  the  Initial 
conditions.  However,  In  the  special  case  of  this  problem,  conditions 
(3.8c,d)  are  consistent  with  (3.9). 

Since  S*"  is  specified,  the  heat  flux  q"*"  must  be  determined  from 
equations  (3.3).  Thus,  using  (3.9)  In  (3.3a)  we  deduce  that 


*  u2 

Po  ch 


(3.10) 


Substituting  (3.10)  Into  (3.3b),  we  have 


*  .2c4  +  x-\».„.  - 

p  ch  [ - z — )  (|»  +  K  (Ji  -  -  q 

°  4x 


(3.11) 


Now,  solving  (3.11)  subject  to  the  Initial  condition  (3.8d),  we  may 
write  the  normalized  solution  In  the  form 


K(e  -  9.)  ,  i  2 

- ^  .  i  [1  -  .TP  (-  (-^^)  Xl]  . 

hq  4  +  11 


[1  -  exp  I-  xl]  . 

q  4+11 


4  ■  [1  -  (-^)  (-  (-^)  xl]  . 


(3.12a) 


(3.12b) 


(3.12c) 


4  +  11 


where  x  la  defined  by  (3.6). 


4  +  It 


To  compare  the  Cosserat  solution  (3.12)  with  the  exact  solution 
recorded  In  [4,  p.  113],  we  rewrite  the  exact  solution  In  the  form 


K(0  -  e^) 


(i  -  9^)  -  Sy  2  sin  x  (i  -  9^)] 

Ti  n-0  (2  n  +  1)“^ 


X  exp  [-  ^  -cj  . 


(3.13) 


Let  US  define  the  average  temperature  9  and  average  temperature 
*  ®vg 

gradient  the  plate  by  the  equations 


*  1  *  3 

9  -  9  -  r  /  (9  -  9  )  de  , 

avg  o  h_i/2 

*  h/2  ^  1  X 

♦avg  “  ^  ®  • 

h^  -h/2  ° 


(3.14a) 


(3.14b) 


Then,  substituting  (3.13)  Into  (3.14)  and  performing  the  Integration,  we 


deduce  the  results 


K(9  -  9  ) 

avg  o 


m 

E 


n-0  (2  n  +  1)^ 


(3.15a) 


-  K 


o 


1  -  ^  ^  [ 
11  n*0 


4  -  (2  n  +  1)  n  (-  1] 


(2  n  +  1) 


-]  exp  [- 


2  n  +  1)  X 


(3.15b) 


Because  the  quantities  6  and  (f)  In  the  Cosserat  solution  correspond 

4t  ig 

to  @  and  (b  .we  have  plotted  each  of  these  In  Figure  B.3.  The 
avg  ^avg 

solid  lines  correspond  to  normalized  values  of  6  and  (|>  and  the  dashed 

it  it 

lines  correspond  to  normalized  values  of  6  and  6  .  The  comparison 

avg  ^avg 

for  all  values  of  x  seems  quite  acceptable. 


Problem  3;  For  this  problem,  we  consider  a  plate  of  thickness  2h.  The 
heat  flux  Is  specified  appropriately  for  radiation  from  both  the  top  and 
bottom  surfaces  and  the  plate  Is  Initially  at  a  uniform  temperature 


9q  +  V.  Mathematically,  these  conditions  are  characterized  by 


q"^  -  KH(0'^  -  0^)  ,  q~  -  -  KH(0"  -  0^)  , 


(3.16a,b) 


0  •  0  +  V  ,  and  ^  «  0  at  t  “  0  , 

o 


(3.16c,d) 


where  H  Is  a  constant  specifying  thermal  radiation  from  the  major 
surfaces.  First,  we  will  solve  the  problem  as  It  Is  formulated  In 
(3.16)  and  second  we  will  obtain  a  more  accurate  solution  by  exploiting 
the  symmetry  about  the  center  plane. 


For  the  first  solution,  we  substitute  (2.161,j)  and  (3.16a,b)  into 
equations  (3.3)  and  then  replace  h  by  2h  to  obtain 


j  ch  § 


-  KH(0  -  0^)  , 


*  2 
4  p  ch 


0  -  -  K(1  +  Hh)  0 


(3.17a,b) 


Using  the  Initial  conditions  (3.16c,d),  the  solution  of  (3.17)  becomes 


0-0 


exp  (-  Hh  x)  ,  0-0  , 


(3.18a,b) 
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FIGURE  D.3  VALUES  OF  THE  NORMALIZED  AVERAGE  TEMPERATURE 
(K(fl  -  AND  AVERAGE  TEMPERATURE  GRADIENT 

[-K  0/q‘]  FOR  A  PLATE  OF  THICKNESS  h  WITH  HEAT  FLUX 
q'  (ENTERING  THE  PLATE)  AT  THE  BOTTOM  SURFACE 
AND  THE  TEMPERATURE  6*  =  0^  SPECIFIED  ON  THE  TOP 
SURFACE 

Initially,  the  temperature  9~  ^  6^  at  the  bottom  surface.  The 
dashed  lines  ate  the  exact  solution  and  the  solid  lines  are  the 
Cosserat  solution.  T  =  Kt/p^  ch  . 
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< 


s  % 


where  t  Is  again  defined  by  (3.6).  To  compare  the  Cosserat  solution 
(3.1R)  with  the  exact  solution  recorded  In  [4,  p.  122],  we  rewrite  the 
exact  solution  In  the  form 


0-0  «  2  Hh  cos  ( — r — )  sec  a  - 

o„  ^h'  n  ,2. 

- - - -  2  - - - - -  exp  (-  t) 

^  n-1  [(Hh)^'  +  Rh  +  a^] 


(3.19a) 


a  tan  a 
n  n 


(3.19b) 


where  are  the  positive  roots  of  equation  (3.19b)  and  where  0*^  ••  0 
locates  the  center  of  the  plate,  0  <■  h  locates  the  top  surface,  and 

0^  ■  -  h  locates  the  bottom  surface.  Replacing  h  by  2h  In  (3.14)  and 
using  (3.19),  we  deduce  the  expressions 


cm 


0  -  0 
avg  o 

V 


s  ^ - 57  exp  (-  a  t) 

n-1  a  [(Hh)'^  +  Hh  +  o'^]  " 


(3.20a) 


(3.20b) 


Comparing  (3.18b)  with  (3.20b),  we  see  that  the  Cosserat  theory  predicts 
the  correct  value  for  the  average  temperature  gradient.  To  compare  the 
prediction  of  the  average  temperature,  we  have  plotted  (3.18a)  as  the 
solid  lines  and  (3.20a)  as  the  dashed  lines  In  Figure  D.4,  for  three 
values  of  the  normalized  radiation  coefficient  Hh.  From  Figure  0.4,  we 
observe  that  for  small  values  of  Hh,  the  Cosserat  theory  predicts 
accurate  results  whereas  for  large  values  of  Hh,  Its  does  not.  This  Is 
because  for  small  values  of  Hh,  heat  Is  radiated  slowly  away  from  the 
major  surfaces  of  the  plate,  so  that  the  temperature  through  the  thick¬ 
ness  of  the  plate  Is  nearly  uniform,  as  predicted  by  (3.18b).  However, 
for  large  values  of  Hh,  heat  is  radiated  rapidly  away  from  the  plate  and 
the  through-the-thlckness  temperature  gradient  may  be  steep. 

Mathematically,  we  may  exploit  the  symmetry  In  the  problem  stated 
above  and  thus  confine  attention  only  to  the  upper  half  of  the  plate. 


D-20 


Therefore,  for  this  second  solution  we  consider  a  plate  of  thickness  h. 
The  heat  flux  Is  specified  appropriately  for  radiation  from  the  top 
surface,  the  bottom  surface  (which  corresponds  to  the  center  surface  of 
the  plate  of  thickness  2h)  Is  Insulated,  and  the  plate  Is  Initially  at  a 
uniform  temperature  6^  +  V.  These  conditions  are  characterized  by 
(3.16)  with  (3.16b)  replaced  by 


(3.21) 


Ac  this  point.  It  Is  Important  to  note  that  although  the  exact  solutions 
of  the  two  problems  considered  here  are  Identical,  the  Cosserat  solution 
of  Che  second  problem  will  be  more  accurate  than  the  Cosserat  solution 
of  Che  first  problem.  This  la  because  Che  solution  of  the  second 
problem  admits  a  nonzero  temperature  gradient  In  the  top  half  of  the 
plate. 


Substituting  (2.161),  (3.16a),  and  (3.21)  Into  equations  (3.3),  we 
obtain 


Po  ch  0  - 


-  1CH(0  +  2  » 


(3.22a) 


*  ^2 
ch 


-  -  KH  (0  +  j  -  0  )  -  K 


(3.22b) 


In  their  present  form  these  equations  are  coupled.  However,  by  solving 
(3.22a)  for  ^  and  substituting  the  result  Into  (3.22b),  we  may  define 


0-0 


-  f(x) 


(3.23) 


and  write 


h_i 

■  V  " 


c/-  ^  2  df 

Hh  dx 


dx 


(3.24a,b) 
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where  B  and  C  are  constants  defined  by 


2  It*  2 

B-ii  +Hh+-^Hh  ,  C-it  Hh  , 


(3.25a,b) 


and  T  Is  defined  by  (3.6).  Using  (3.6),  (3.23),  and  (3.24a),  the 
Initial  conditions  (3.16c,d)  become 


f  -  1 

^  ^  ’  dx 


-  Hh  ,  at  T  *  0 


(3.26a,b) 


Solving  (3.24b)  subject  to  the  conditions  (3.26),  the  Cosserat  solution 
may  be  written  In  the  form 


0-9 


_  h  (fc  _ 


exp  (-  OjT)  +  exp  (- 


(3.27a) 


^  -  2  [A^  exp  (-  o^t)  +  A^  exp  (-  a^x)] 


-  exp  (-  Oj^x)  <32  ^2*^^^  *  (3.27b) 

where  the  constants  A^,  A2,  02  are  given  by 


02  -  Hh  Hh  -  Oj^ 

^1  "  CT2  -  *  ^2  "  <^2  -  o^  ’ 


Oj  -  Y  [B  -  (B^  -  4  ,  02 


(3.28a,b) 


1  2  1/2 
Y  [B  +  (B  -  4  C)"'^] 


(3.28,c,d) 


Replacing  0^  In  (3.19a)  by  h/2  +  0^,  we  may  write  the  exact 
solution  for  the  top  half  of  the  plate  as 


2  Hh  cos  j 

i  a  (h  +  2  9^)  1 
)— - ;-r - 1 

[  sec  a 

1 

[  2  h  J 

’  n 

[(Hh)^  +  Hh  +  a  ] 


exp  (-  x)  , 


(3.29) 


IVIVI' 


where  are  the  positive  roots  of  (3.19b),  and  where  9^  ■  h/2  locates 
the  top  surface  and  6  >  -  h/2  locates  the  bottom  surface  (which 

corresponds  to  the  center  surface  of  the  plate  of  thickness  2h) . 

Substituting  (3.29)  into  the  definitions  (3.14a,b),  we  obtain  the  result 

* 

(3.20a)  for  the  average  temperature  ^nd  the  result 


«  Hh  fa  sin  a  +  2  (cos  a  -  1)1  sec  a  _ 

n  n  n  n 

12  2  - T - 5 - 5 -  exp  (-  a  t)  . 

n-1  a  [(Hh)'^  +  Hh  +  a  ] 

"  “  (3.30) 


Figure  D.5  compares  values  of  the  average  temperature  with  (3.27a) 
plotted  as  the  solid  lines  and  (3.20a)  plotted  as  the  dashed  lines. 
Similarly,  Figure  D.6  compares  values  of  the  average  temperature 
gradient  with  (3.27b)  plotted  as  the  solid  lines  and  (3.30)  plotted  as 
the  dashed  lines.  From  Figures  D.4  and  D.5,  we  observe  that  modeling 
only  the  upper  half  of  the  plate  produces  a  significant  improvement  in 
the  prediction  of  the  average  temperature  for  the  higher  values  of  Hh. 
Also,  we  observe  from  Figures  0.5  and  0.6  that  for  Hh  >  0.1,  the 
Cosserat  and  exact  solutions  are  nearly  identical  and  the  average 
temperature  gradient  remains  relatively  small. 


Circular  Cylindrical  Shells 


In  this  section,  we  investigate  the  validity  of  the  Cosserat  theory 

in  the  thick-shell  limit  by  considering  heat  conduction  in  a  circular 

cylindrical  shell  and  taking  the  limit  of  a  solid  cylinder.  Specifically, 

consider  a  circular  cylindrical  shell  of  uniform  thickness  h  and  mean 

radius  R.  Let  e.(i  ~  1,2,3)  be  a  set  of  fixed  Cartesian  base  vectors 
~1 

and  let  e'  be  base  vectors  of  a  polar  coordinate  system  with  polar  angle 

* 

Y  defined  by 


Although  this  coordinate  system  is  unconventional,  it  is  chosen  because 
it  yields  convenient  relations  between  and 
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FIGURE  D.6  VALUES  OF  THE  NORMALIZED  AVERAGE  TEMPERATURE  GRADIENT 
(-h  0/V]  FOR  A  PLATE  OF  THICKNESS  h  WITH  RADIATION  FROM 
THE  TOP  SURFACE  (i.e.,  q'^  =  KH  W"  -  0^)] ,  ZERO  HEAT  FLUX  ON  THE 
BOTTOM  SURFACE,  AND  UNIFORM  INITIAL  TEMPERATURE  «  =  +  V 

The  dashed  lines  are  the  exact  solution  and  the  solid  lines  are  the  Cosserat  solution. 
T  =  Kt/p  *  ch^ . 
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15^0 


(4.6) 


,2  9  .  sLe  + _ 

»  o  9  2  2 

a  R  a  Y 


Here,  we  consider  the  problem  for  which  the  heat  flux  on  the  outer 
surface  Is  constant,  the  Inner  surface  Is  Insulated,  external  entropy 
supply  Is  neglected,  and  the  shell  Is  Initially  at  uniform  temperature 
6^.  Hence,  the  conditions  (3.1)  and  (3.4)  hold  and  equations  (4.5) 
reduce  to 


Pq  ch  0  -  -  (l  +  1^]  q"*"  ,  (4.7a) 

p*  ch  ;  -  -  K  ^  -  Y  (1  +  ^)  q"^  .  (4.7b) 

Integrating  (4.7)  subject  to  the  Initial  conditions  (3.4c,d),  we  obtain 


K(0  -  0^)  . 

- ^  .  (1  +  ^  ,  (4.8a) 

h  q 

-  ^  -  Y  (1  +  -  exp  (-  ^  t)!  .  (4.8b) 

H 

where  t  Is  defined  by  (3.6).  Notice  that  In  the  thln-shell  limit 
(R/h  -*■  •),  the  solution  (4.8)  approaches  the  plate  solution  (3.5).  In 
the  thick-shell  limit  of  a  solid  cylinder  for  which*  R  ■  h/2,  the  right- 
hand  side  of  (4.8a)  becomes  2t,  which  Is  consistent  with  the  exact 
solution  [4,  p.  203].  Using  (2.15b),  the  Cosserat  solution  (4.8)  with 
R  ■  h/2.  Is  plotted  In  Figure  D.7  together  with  the  exact  solution  for 
various  values  of  the  time  parameter  t.  The  dashed  lines  In  Figure  D.7 
have  been  taken  directly  from  [4,  Figure  25,  p.  203]  and  r  Is  the  radial 


From  [2],  we  recall  that  generally  Is  required  to  be  positive. 

Although  the  quantity  vanishes  when  0^  •  -  h/2  and  R  ■  h/2,  this 

poses  no  particular  difficulty  in  the  problems  considered  here. 
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FIGURE  D.7  NORMALIZED  TEMPERATURE  IN  A  SOLID  CIRCULAR 
CYLINDER  OF  RADIUS  h,  WITH  CONSTANT  HEAT 
FLUX  q*  (OUT  OF  THE  CYLINDER)  AT  THE  SURFACE, 
AND  UNIFORM  INITIAL  TEMPERATURE  6=0^ 

The  numbers  on  the  curved  lines  are  values  of  t  =  Kt/p^  ch^ . 
The  dashed  lines  are  the  exact  solution  and  the  solid  lines  are 
the  Cosserat  solution. 
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coordinate,  with  r  >  0  locating  the  center  of  the  cylinder  and  r  -  h 
locating  the  outer  surface. 

In  view  of  the  form  of  solution  (4.8),  It  Is  obvious  that  the  long¬ 
time  temperature  Is  dominated  by  the  term  (4.8a).  The  fact  that  the 
coefficient  of  t  In  (4.8a)  yields  the  correct  result  even  In  the  thick- 
shell  limit  suggests  that  the  Cosserat  theory  retains  the  Important 
geometrical  features  of  the  shell.  In  this  regard.  It  Is  worth  men¬ 
tioning  that  the  result  (4.8a)  could  be  obtained  using  an  engineering 
approach  In  which  the  temperature  In  the  shell  Is  assumed  to  be  uniform 
and  the  energy  entering  the  outer  surface  Is  equated  with  the  Increase 
In  Internal  energy.  It  Is  also  worth  mentioning  that  the  solution  of 
the  more  accurate  equation  (15)  of  [1]  yields  a  long-time  solution  of 
the  form 


K(9*  -  9^) 


Fslnh  +  cosh 

L  ^  (fe)  J" 


(4.9) 


In  the  thln-shell  limit  (4.9)  yields  the  correct  result,  but  In  the 
thick-shell  limit  it  yields  the  result  (2.313  t),  which  Is  Incorrect. 
Thus,  even  though  equation  (15)  of  [1]  Is  more  complicated  than 
equations  (4.8),  It  does  not  necessarily  produce  a  better  result. 

To  further  examine  the  validity  of  the  constitutive  equations 
(2.14b,c)  and  the  specifications  (2.18a,c),  we  consider  a  simple  problem 
for  which  the  Laplaclan  operators  In  (4.5)  do  not  vanish.  Specifically, 
consider  the  steady-state  problem  of  uniform  heat  conduction  in  the 
constant  £3  direction  for  which  the  three-dimensional  solution  is  given 
by 


S  “  Q  £3  »  ®  ■  (Qq  ■  ^  cos  y)  ■  cos  y)  (4.10a,b) 

* 

where  ^  Is  the  three-dimensional  heat  conduction  vector  and  Q  Is  a 
constant.  Using  (4.1)  and  (4.10),  we  realize  that 
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Consequently,  in  the  absence  of  external  entropy  supply,  the  steady- 
state  solution  of  (4.5)  becomes 

2  -1 

0  ■  9  -  ^  cos  y  ,  (>--(!+ - r)  ^  cos  y  .  (4.12a,b) 

°  12TC  ^ 

Now,  with  the  help  of  (2.15b)  we  may  compare  the  exact  result  (4.10b) 
with  the  Cosserat  result  (4.12)  to  conclude  that  the  average  temperature 
is  predicted  exactly.  Further,  the  prediction  of  the  average  tempera¬ 
ture  gradient  ^  la  very  accurate  in  the  thln-shell  limit  (R/h  <»)  and 
is  only  25%  low  in  the  thick-shell  limit  (R/h  -*•  1/2). 

5.  Spherical  Shells 

The  spherical  shell  geometry  is  considered  here  mainly  because  it 
is  one  of  the  simplest  geometries  in  which  it  is  possible  to  investigate 
the  differences  between  the  constitutive  assumptions  (2.12)  and  (2.14). 
Three  problems  of  a  spherical  shell  of  constant  thickness  h  and  radius  R 
are  considered.  For  the  first  two  problems,  we  consider  the  thick-shell 
limit  of  a  solid  sphere  and  discuss  the  differences  between  assumptions 
(2.12)  and  (2.14).  For  the  third  problem,  we  consider  the  transition 
from  a  thick  shell  to  a  thin  shell. 

With  reference  to  the  Cartesian  base  vectors  Introduced  in 

Section  4,  we  let  e' '  be  base  vectors  of  a  spherical  coordinate  system 
~1 

with  polar  angle  y  (0  <  y  <  2x)  measured  from  the  e  -  e  plane  and 

polar  angle  a  (-  —  <  a<-:r]  measured  from  the  e  -  e  plane  such  that 
'•  2  •  2  '^1  '^2 

e"  •  -  sin  Y  e  +  cos  y  e  ,  (5.1a) 

~1  ~2 

*Although  this  coordinate  system  is  unconventional,  it  is  chosen  because 
it  yields  convenient  relations  between  and  fej'* 


(5.1b) 


e''  ■  -  sin  a  (cos  v  e,  +  sin  v  e.)  +  cos  a  e*  , 
^2  *^1  ^^3 


~3 


cos  a  (cos  Y  Y  £.2^  °  ~3 


(5.1c) 


Now,  points  on  the  reference  surface  of  the  shell  may  be  located  by 
the  position  vector  g  given  by 


R  •  R  el'  ,  0 

~  ~3 


(5.2a,b,c) 


i?:-: 


1  o 

where  we  have  Identified  the  coordinates  Q‘- ,  0  with  y  and  <j, 
respectively.  Using  the  definitions  In  [2]  and  In  Section  2,  the 
relevant  geometrical  properties  of  the  spherical  surface  may  be  recorded 
as 


„2  11 

A  ■  R  cos  o  ,  A  ■ 


n2  2 

R  cos  0 


«1  -  ®2 


.12  ^  22  1 

,  A  -  0  ,  A  "  , 

R 

(5.3a,b,c,d) 


,  all  other  B*  -  0  , 

R  p 


(5.3e,f) 


r.1  r.1 

^12  “  ^21 


,2 

11 


tan  0  ,  r,i  *  sin  o  cos  a 


all  other  F  _  •  0 

ap 


(5.3g,h,i) 
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<jl 

'hi 

r'-l 

m 

0] 


IjfJ 


m 

*s 


^  (1  +  »  )  ,  .  Po  9o  -1  -  mi  +  q  -  (f)(l  -  q 

t  12  R 

2  3 

-  Kh  (1  +  ^  7^  (>  ,  (5.5b) 

12  R*^ 


where  the  Laplaclan  operator  7  0  for  the  spherical  geometry  Is  given  by 


2  1  9  1_  a^  9  tan  a  W 

®  “  _2  2  .  2  „2  .  2  ~  _2  acr  ■ 

R  cos  a  a  Y  R  a  o  R 


(5.6) 


For  the  first  problem,  the  heat  flux  on  the  outer  surface  Is 
constant,  the  Inner  surface  Is  Insulated,  external  entropy  supply  Is 
neglected,  and  the  shell  Is  Initially  at  uniform  temperature  9^.  Hence, 
the  conditions  (3.1)  and  (3.4)  hold  and  equations  (5.5)  reduce  to 


p*  ch  1  +  §  -  -  (1  +  q+  . 

°  12  R^ 


(5.7a) 


*  2 

ch  u2  .  1  u  t  M  k2 

-2_-  (1  +-S!— r)  ♦  -  -  i  (1  q+  -  K  (1  +-Si— r)  ♦  .  (5.7b) 


Integrating  (5.7)  subject  to  the  Initial  conditions  (3.4c,d),  we  obtain 


Kca  -  9„)  [1  . 


(5.8a) 


K  4 
+ 
q 


(1  + 


2R^  „  ,  2  .. 

— 2 -  [1  -  exp  (-  X  t)I  , 


(5.8b) 


2(1+-^) 


where  x  Is  defined  by  (3.6).  Notice  that  In  the  thln-shell  limit 
(R/h  «)  the  solution  (5.8)  approaches  the  plate  solution  (3.5).  In 
the  thick-shell  limit  of  a  solid  sphere  for  which  R  ■  h/2,  the  rlght- 


►*  *1 

lit 

y 


hand  side  of  (5. 8a)  becomes  3t»  which  Is  consistent  with  the  exact 
solution  p.  242].  Using  (2.15b),  the  Cosserat  solution  (5.8)  with 
R  -  h/2  Is  plotted  In  Figure  D.8  together  with  the  exact  solution  for 
various  values  of  the  time  parameter  t.  The  dashed  lines  In  Figure  D.8 
have  been  taken  directly  from  [4,  Figure  31,  p.  242]  and  r  is  the  radial 
coordinate  with  r  -  0  locating  the  center  of  the  sphere  and  r  >  h 
locating  the  outer  surface. 

From  Figure  D.8,  we  observe  that  for  long  time  periods  the  value  of 
the  average  temperature  gradient  predicted  by  the  Cosserat  theory  Is 
substantially  larger  than  the  exact  value.  However,  this  Is  not 
particularly  Important  because  for  long  times  the  temperature  Is 
dominated  by  the  term  (5.8a).  To  exhibit  this,  we  have  used  (2.161,j) 
together  with  (5.8)  to  plot  In  Figure  D.9  the  temperature  on  the  outer 
surface  and  at  the  center  of  the  solid  sphere.  The  dashed  lines  In 
Figure  D.9  represent  the  exact  solution  [4,  p.  242].  For  short  times, 
the  Cosserat  theory  predicts  the  Incorrect  result  that  the  center 
temperature  of  the  sphere  drops.  This  Is  a  consequence  of  the  over- 
prediction  of  the  average  temperature  gradient.  For  long  times,  the 
lines  In  Figure  D.9  are  parallel  and  the  relative  error  diminishes  to 
zero.  Thlq  is  because  the  prediction  (5.8a)  Is  exact  In  the  thick-shell 
limit.  In  this  regard.  It  Is  worth  mentioning  that  the  result  (5.8a) 
could  be  obtained  using  the  engineering  approach  described  In  Section  4. 
It  Is  also  worth  mentioning  that  the  more  accurate  equation  (15)  of  [1] 
yields  a  long-time  solution  of  the  form* 


K(9*  -  9^)  "slnh  +  cosh  (^)' 
h  q"^  _  slnh  (^) 


(5.9) 


*The  solution  In  Appendix  A  of  [1]  should  be  written  In  a  form  that  has 
a  linear  term  In  time. 
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FIGURE  D.8  NORMALIZED  TEMPERATURE  IN  A  SOLID  SPHERE 
OF  RADIUS  h,  WITH  CONSTANT  HEAT  FLUX  q"" 
(OUT  OF  THE  SPHERE)  AT  THE  SURFACE,  AND 
UNIFORM  INITIAL  TEMPERATURE  fl 

The  numbers  on  the  curves  are  values  of  t  =  Kt/p*  ch^ . 
The  dashed  lines  are  the  exact  solution  and  the  solid  lines 
are  the  Cosserat  solution. 
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FIGURE  0.9 


NORMALIZED  TEMPERATURE  I-KIfl*  -  0^)/hq*)  ATTHE  SURFACE 
AND  *  -  O^lW]  AT  THE  CENTER  OF  A  SOLID  SPHERE  OF 


RADIUS  h,  WITH  CONSTANT  HEAT  FLUX  q*  (OUT  OF  THE  SPHERE) 
AT  THE  SURFACE,  AND  UNIFORM  INITIAL  TEMPERATURE  ^ 


The  dashed  lines  are  the  exact  solution  and  the  solid  lines  are  the  Cosserat 


solution.  T  »  Kt/o  ch^ 


D-36 


*- w"- . 


■f  / 


In  the  chin-shell  limit  (5 >9)  yields  the  correct  result,  but  In  the 
thick-shell  limit  It  yields  Che  result  (4.075  x),  which  Is  Incorrect 


He  are  now  In  a  position  to  comment  on  the  differences  between  the 
constitutive  equations  (2.12)  and  (2.14).  If  (2.12a)  were  used  Instead 
of  (2.14a),  then  the  average  temperature  would  be  given  by 

K(0  -  0„)  .  , 

-  -  (1  +  If)'  ^  (S-IO) 

h  q 

Instead  of  (5.8a).  This  would  yield  the  Incorrect  result  4t  In  Che 
chick-shell  limit.  Similarly,  If  (2.12d,e)  were  used  Instead  of 
(2.14d,e),  then  the  long  time  value  of  would  be 

q 

which  produces  a  larger  error  than  that  associated  with  (5.8b)  In  the 
chick-shell  limit. 

To  further  examine  the  validity  of  constitutive  equations  (2.14b,c) 
and  Che  specifications  (2.18a,c),  we  consider  the  steady-state  problem 
of  uniform  heat  conduction  In  the  constant  direction  for  which  the 
three-dimensional  solution  Is  given  by 

£  *  0  "^®o”^  ^  ^  ’ 

(5.12a,b) 

* 

where  £  Is  the  three-dimensional  heat  conduction  vector  and  0  is  a 
constant.  TJslng  (5.1)  and  (5.12),  we  realize  that 

q^-q  ■£  ‘ej'-Q  cos  y  cos  a  .  (5.13) 

Consequently,  In  Che  absence  of  external  entropy  supply  the  steady-state 
solution  of  (5.5)  becomes 


I)-37 


cos  Y 


(5.14a,b) 


©  -  0 

®  ®o  K 


,  ♦  ■  "  K  y  ®  » 


which  Is  an  exact  result  valid  for  both  the  thin-  and  thick-shell 
limits. 

Finally,  to  examine  the  transition  from  a  thin-shell  to  a  thick- 
shell,  we  consider  the  steady-state  problem  where  the  temperature  on 
the  outer  surface  is  specified  to  be  the  constant  value  ©^  and  the  heat 
flux  q~  on  the  inner  surface  is  constant.  Thus,  using  (2.161)  we 
require 


©  +  Y4*©^  ,  q  -  constant  . 


(5.15a,b) 


In  the  absence  of  external  entropy  supply,  the  steady-state  solution  of 
(5.S)  becomes 


K(©-©o)  (1-^)' 


2(1  + 


fl  - 

(i  +  -*W) 


(5.16a,b) 


(5.16c) 


It  can  be  shown  that  the  exact  solution  [4,  p.  247]  may  be  written  in 
the  form 


ic(e*-9„)  Ci-W)^(i-V-) 


(5.17) 


2(1  +  ^)(1  .  f-) 


and  that  (5.16c)  is  an  exact  result.  Now  to  compare  the  predictions 
(5.16)  with  the  exact  solution  (5.17),  we  have  used  (2.15b)  to  plot 
(5.16)  as  the  solid  lines  in  Figure  D.IO  and  have  used  (5.17)  to  plot 
the  dashed  lines  in  Figure  D.IO  for  three  values  of  R/h.  The  results  in 
Figure  0.10  show  again  that  the  Cosserat  predictions  are  good  even  for  a 
fairly  thick  shell  (R/h  ■  1). 
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FIGURE  D.10  THE  NORMALIZED  STEADY  STATE  TEMPERATURE  IN 
A  SPHERICAL  SHELL  WITH  CONSTANT  THICKNESS  h 
AND  MEAN  RADIUS  R 

The  heat  flux  c|  (entering  the  shell)  is  specified  on  the  inner  surface 
and  the  temperature  $*  *  0^  is  specified  on  the  outer  surface.  The 
dashed  lines  are  the  exact  solution  and  the  solid  lines  are  the  Cosserat 
solution. 
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Conical  Shell 


In  the  previous  sections,  we  have  solved  a  number  of  problems  for 
plates,  circular  cylindrical  shells,  and  spherical  shells  to  develop 
confidence  that  the  Cosserat  theory  can  predict  relatively  accurate 
results  for  both  the  thln-shell  limit  (which  models  the  base  of  a 
conical  shell)  and  the  thlck'-shell  limit  (which  models  the  tip  of  a 
conical  shell).  Here,  we  confine  attention  to  a  conical  shell  with 
constant  thickness  h  and  locate  points  on  the  conical  surface  by 


R  -  p  R  e  ’  +  R  e^  ,  0^  -  R  ,  9^  -  y  , 


(6.1a,b,c) 


where  R  Is  the  radial  coordinate,  y  Is  the  polar  angle,  e'  are  defined 

~1 

by  (4.1),  p  Is  a  constant  related  to  the  cone  angle  (see  Figure  D.l), 
and  we  have  Identified  the  coordinates  9^  and  9^  with  R  and  y, 
respectively.  Using  the  definitions  In  [2]  and  In  Section  2,  the 
relevant  geometrical  properties  of  the  conical  surface  may  be  recorded 


.1/2  .  .2^/2  .11  1  12  ..  .22  1 
A  -  R(1  +  P  )  ,  A  »  - ST*  .A  -  0  ,  A  -  “T  , 

(1  +  P'^)  R'^ 


(6.2a,b,c,d) 


2  6  (Z 

B,  - o'  T  H  »  «11  other  B®  •  0  , 

^  R(1  +  ° 


(6.2e,f) 


R  _2  «2  1  . .  .  _a  .. 

“■“2^  •  ^12  -  ^21  "  R  '  ^ap  “  °  • 


(6.2g,h) 


Substituting  (6.2)  Into  (2.16),  we  have 


*  ^  +  l/2r,  h  p 

po "  po  ^  ®  ^  - TTn 

2R(1  +  p  )  ' 


1  »■  *1/2,  h  p 

]  »  ®  “  A  [1  -  2  1/2I 

2R(1  +  p  )  ' 


(6.3a,b,c) 


vV 


It  follows  that  the  thermal  equations  (2.19)  become 


I*; 

iM 


I 


*  * 

p  ch  9 


p*  h  9  s  -  [  1  + - — -y  I'/V 

°  °  2R(1  + 


+  K  h  V  9  , 


]  +  [1  - 


n  p  T  - 

-  2  1/2  J 

2R(1  +  p  )  ' 


(6 .4a) 


* 

*^0  I  .  *  ft  :  i  ri  4.  h  S  1  +  i  r,  h  9  1  - 

2  ^  *^o  ®o  ®1  ”  2  2  1/2-1  “2  1-^  2  1/2 J 

n  o  o  1  z  2r(1  +  2R(1  +  p^)^' 


-  K  ♦  +  §-  ♦  , 


(6.4b) 


where  the  Laplacian  operator  7^  0  for  the  conical  geometry  Is  given  by 


2  2 

ri  9i 


(6.5) 


Here,  we  consider  the  problem  for  which  the  heat  flux  on  the  outer 
surface  Is  constant,  all  other  surfaces  are  Insulated,  external  entropy 
supply  Is  neglected,  and  the  shell  Is  initially  at  uniform  temperature 
9q.  Hence,  the  conditions  (3.4)  hold  in  addition  to  the  conditions 


9  -  9(R,t)  ,  a  -  ♦(R.t)  , 


0  ^ 
911  ■  9R 


0  at  R 


(6.6a,b) 


(6.6c,d) 


where  R^  and  R2  are  the  tip  radius  and  base  radius  of  the  shell, 
respectively  (see  Figure  D.l).  Under  these  conditions,  equations  (6.4) 


reduce  to 


Po  ch  9 


ri  ^  h  p  1  ^  Kh  a  ,  a99 

2R(1  +  p^)^^^  ^  R(1  +  p^) 


(6.7a) 


♦ 
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*  u2 

p  ch 
- 5 -  9 


-^1  + 


JlA. 
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To  analyze  these  equations.  It  is  convenient  to  Introduce  the  nondlmen- 
slonal  parameters 
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and  rewrite  them  In  the  form 
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Similarly,  the  Initial  conditions  (3.4c,d)  and  boundary  conditions 
(6.6c,d)  become 
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where  z^  and  Z2  are  the  values  of  z  when  R  equals  R^  and  R2, 
respectively.  At  this  point.  It  Is  of  Interest  to  note  that  In  the 
limit  of  large  P(P  -»■  »),  equations  (6.9)  reduce  to  a  nondlmenslonal  form 
of  (4.7)  for  a  circular  cylindrical  shell,  and  In  the  limit  of  small 
p(p  0),  equations  (6.9)  characterize  a  circular  plate. 

Using  standard  techniques,  the  solution  of  (6.9)  may  be  written  in 
the  form 
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where  fm(z)  are  eigenfunctions  characterized  by 
,  df 

—  -j—  (z  -  a  f  (no  8\im  on  m)  (6.12a) 
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and  where  are  the  nonzero  eigenvalues.  Since  equation  (6.12a)  can 
easily  be  recognized  as  Bessel's  equation  of  order  zero,  the  solution, 
subject  to  the  boundary  conditions  (6.12b),  Is  well  characterized  and 
may  be  written  In  the  form 
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(6.13) 


where  and  are  Bessel  functions  of  the  first  and  second  kind, 
respectively,  of  order  n  and  where  are  the  positive  roots  of  the 
characteristic  equation 

J  (a  z  )Y  (a  z  )  -  Y  (a  z  )J  (a  z  )  -  0  .  (6.14) 
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Further,  the  eigenfunctions  f^,  satisfy  the  orthogonality  conditions 
*2  *2 

/  zfdz*0,  /  zffdz"0  for  (m  ^  n) 

•'  m  ’  m  n 


(6.15a,b) 


Substituting  (6.11)  into  (6.9),  multiplying  the  result  by  z  and 
integrating,  multiplying  the  result  by  zf^  and  integrating,  and  using 
the  orthogonality  condition  (6.15)  and  the  initial  conditions  (6.10a,b), 
we  conclude  that 
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For  later  reference,  we  observe  that  if  the  dependence  on  z  is  neglected 
in  (6.11),  then  (6.11)  has  the  same  form  as  the  solution  (4.8).  This 
means  that  %ie  would  be  essentially  modeling  the  conical  shell  as  an 
"equivalent"  circular  cylindrical  shell  with  "mean"  radius  R/h  *  (z^  + 
22)(1  +  p2)^^2/23.  By  considering  a  specific  example,  it  will  be  shown 
that  making  this  kind  of  engineering  approximation  Introduces  signif¬ 
icant  errors  at  the  tip  of  the  conical  shell. 

Consider  the  specific  example  of  the  conical  shell  drawn  in  Figure 
T).l,  which  has  a  solid  tip.  For  this  shell  we  specify 
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The  minimum  value  of  z  given  by  (6.17a)  Is  specified  by  requiring  Che 
Inner  surface  of  the  shell  to  Just  make  contact  at  the  shell's  dp. 

Using  (6.17)  we  have  solved  for  the  first  twenty  eigenvalues  and  eigen¬ 
functions  and  have  plotted  the  solution  (6.11)  In  Figures  D.ll  and  D.12 
by  normalizing  the  results  by  the  first  terms  In  the  solutions.  Figure 
D.ll  shows  plots  of  ?/Aq(x)  verses  z  for  various  values  of  t,  and 
Figure  D.12  shows  plots  of  ^/Bq(x)  verses  z  for  two  values  of  x.  The 
slight  wavlness  In  these  curves  Is  caused  by  Che  fact  chat  we  have 
approximated  solutions  (6.11a,b)  using  finite  series. 

From  Figure  D.ll,  we  observe  that  for  long  times  the  average 
temperature  Is  relatively  uniform  over  the  shell.  This  Is  because  the 
equivalent-cylinder  solution  Aq(x)  dominates  for  long  times.  However, 
for  short  times  Che  value  of  9  at  the  tip  Is  about  65%  greater  than  that 
predicted  by  Che  equivalent-cylinder  solution.  This  result  can  be 
explained  by  observing  from  (4.8a)  chat  a  thick  cylinder  heats  up  faster 
Chan  a  thin  cylinder.  Thus,  we  would  expect  the  tip  of  the  conical 
shell,  which  Is  thick,  to  heat  up  faster  than  Its  base,  which  Is  thin. 
From  Figure  D.12,  we  observe  that  the  distribution  of  the  average 
temperature  gradient  Is  nearly  constant  with  time.  Also,  the  value  of  $ 
near  the  dp  Is  nearly  65%  greater  than  the  value  predicted  by  the 
equivalent-cylinder  solution. 

To  exhibit  Che  temporal  dependence  of  this  solution  more  clearly, 
we  have  plotted  Aq(x)  and  Bq(x)  In  Figure  D.13.  From  this  figure,  we 
observe  that  Bg(T)  reaches  Its  maximum  value  In  a  relatively  short  time. 
Recalling  from  [2]  Chat  Che  average  temperature  gradient  Is  related  to 
the  thermal  bending  moment  In  the  shell,  this  means  that  the  bending  In 
a  conical  shell  uader  this  load  will  be  quite  severe  at  its  dp  and  Che 
full  effect  of  the  load  will  be  felt  In  a  relatively  short  time. 
Consequently,  Che  dp  of  Che  conical  shell  should  be  particularly 
vulnerable  to  this  type  of  thermal  load. 


FIGURE  D.11  NORMALIZED  AVERAGE  TEMPERATURE  IN  A  CONICAL  SHELL 
OF  THICKNESS  h. 

The  numbers  on  the  curves  are  values  of  t  »  Kt/p*  ch^. 
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In  this  paper,  we  have  focused  attention  on  analyzing  heat 
conduction  In  a  rigid  conical  shell  (Figure  D.l).  The  conical  shell  Is 
particularly  Interesting  because  It  has  a  converging  geometry,  so  that 
the  shell  near  Its  tip  Is  necessarily  "thick”  even  though  the  shell  near 
Its  base  may  be  "thin.”  Further,  the  heat  conduction  equation  Is  not 
separable  for  the  conical  geometry,  and  hence  It  Is  exceedingly 
difficult  to  obtain  exact  solutions.  Ue  have  chosen  to  model  the  shell 
with  the  theory  of  a  Cosserat  surface  to  determine  the  average  (through- 
the-thlckness)  temperature  and  temperature  gradient  In  the  shell 
directly. 

A  number  of  problems  of  plates,  circular  cylindrical  shells,  and 
spherical  shells  are  considered  and  the  solutions  are  compared  with 
exact  solutions  to  develop  confidence  In  the  Cosserat  theory.  Within 
the  context  of  this  theory.  It  Is  usually  assumed  that  constitutive 
equations  for  shells  have  the  same  form  as  those  for  plates.  Here,  It 
Is  shown  that  to  predict  relatively  accurate  results  In  the  thick-shell 
limit.  It  Is  necessary  to  generalize  these  constitutive  equations  to 
Include  certain  geometrical  features  of  the  shell.  The  generalized 
constitutive  equations  are  developed  here  In  a  consistent  manner  and 
tested  In  the  thick-shell  limit.  The  tests  Include  problems  where  the 
temperature  fields  9  and  ^  are  functions  of  time  only  so  that  their 
Laplaclan  vanishes,  as  well  as  problems  where  they  are  functions  of 
space  only  and  their  Laplaclan  does  not  vanish.  In  all  cases, 
satisfactory  results  are  predicted  even  In  the  thick-shell  limit. 

Finally,  a  problem  of  transient  heat  conduction  In  a  conical  shell, 
which  does  not  have  an  exact  solution.  Is  solved  analytically  using  the 
Cosserat  theory.  It  Is  shown  that  both  the  average  temperature  and 
temperature  gradient  have  values  near  the  tip  that  are  about  65Z  greater 
than  those  predicted  by  an  approximate  equivalent-cylinder  solution. 
Also,  It  Is  shown  that  the  thermal  bending  moment  produced  by  the 
average  temperature  gradient  Is  quite  severe  near  the  tip  and  It  attains 
Its  maximum  value  In  a  relatively  short  time. 
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